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ARBITRARY  ASPECT  RATIO,  SUPERCAVITATING 
HYDROFOILS  OPERATING  NEAR  THE 
FREE  WATER  SURFACE 
By  Virgil  E.  Johnson,  Jr. 

SUMMARY 


Theoretical  expressions  for  the  lift,  drag,  and  center  of  pressure 
of  supercavitating  hydrofoils  of  arbitrary  section,  aspect  ratio,  and 
depth  of  submersion  are  developed  for  the  zero-cavitation-number  case. 

In  order  to  verify  and  supplement  the  theoretical  investigation, 
an  experimental  investigation  was  made  in  Langley  tank  no.  2  of  two 
aspect-ratio-1  hydrofoils,  one  with  a  flat  and  one  with  a  cambered  lower 
surface.   Zero  cavitation  number  was  obtained  in  the  tank  by  operating 
the  hydrofoils  near  the  water  surface  so  that  their  upper  surfaces  were 
completely  ventilated.   Data  were  also  obtained  at  cavitation  numbers 
greater  than  zero  for  flows  with  vapor-filled  cavities  and  with  air  arti- 
ficially introduced  to  simulate  the  same  flows  at  lower  speeds.   The  tests 
covered  a  range  of  water  speeds  up  to  80  feet  per  second,  angles  of 
attack  up  to  20°,  and  depths  of  submersion  from  0-  to  85-percent  chord. 

The  theoretical  force  and  moment  coefficients  agreed  with  the  data 
obtained  in  the  zero-cavitation-number  tests  within  an  accuracy  of  about 
3  percent.   The  ventilated  force  and  moment  coefficients  of  both  models 
were  independent  of  speed.   As  predicted  by  the  linear  theory  of  Tulin 
and  Burkart,  the  cambered  hydrofoil  lift-drag  ratio  was  superior  to  that 
for  the  flat  plate  at  the  same  lift  coefficient.   Similar  to  the  influ- 
ence of  camber  on  a  fully  wetted  airfoil,  the  influence  of  camber  in 
ventilated  flow  was  to  effectively  increase  the  angle  of  attack  without 
appreciably  changing  the  lift-curve  slope.   Comparison  of  the  present 
results  with  the  results  of  other  investigations  revealed  that  a  sharp 
leading  edge  is  necessary  for  good  efficiency  in  supercavitating  flow. 
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The  nonzero-cavltation-number  tests  confirmed  the  finding  of  others 
that  the  lift  of  a  hydrofoil  operating  in  a  finite  cavity  depends  only 
on  the  cavitation  nijmber  and  is  independent  of  whether  the  gas  in  the 
cavity  is  air  or  vapor. 


INTRODUCTION 


The  desirability  of  using  aiixiliary  lifting  surfaces  such  as  hydro- 
foils or  hydroskis  for  reducing  seaplane-hull  loads  and  improving  rough- 
water  performance  is  well  established.   Although  the  higher  aspect-ratio 
submerged  hydrofoils  with  t]ieir  accompanying  high  efficiencies  could  be 
the  most  desirable  lifting  surfaces  to  be  used,  only  the  low-aspect-ratio 
planing  hydroski  has  so  far  been  successfully  applied  as  landing  gear  to 
modern  high-speed  aircraft.   The  reason  for  this  is  that  the  conventional 
hydrofoil  presents  problems  not  experienced  by  a  hydroski. 

As  a  hydrofoil  nears  the  free  water  surface  (during  a  take-off  run) 
the  low  pressure  side  of  the  hydrofoil  becomes  ventilated  from  the  atmos- 
phere.  This  phenomenon  results  in  a  severe  and  usually  abrupt  loss  in 
lift  and  reduction  in  the  lift-drag  ratio.   For  conventional  airfoil 
sections,  the  loss  in  lift  may  exceed  75  percent.  The  speed  at  the 
inception  of  ventilation  depends  on  the  angle  of  attack  and  depth  of 
submersion;  but,  except  for  very  small  angles  of  attack  and  relatively 
low  take-off  speeds,  the  speed  at  Inception  is  usually  well  below  the 
take-off  speed  of  the  aircraft. 

Even  if  the  ventilation  problem  is  overcome  by  using  small  angles 
of  attack  and  incorporating  fences  or  other  devices  for  suppressing 
ventilation,  the  onset  of  cavitation  presents  a  second  deterrent  to  the 
use  of  conventional  hydrofoils  at  high  speeds.   The  loss  in  lift  accom- 
panying cavitation  of  conventional  airfoil  sections  is  not  abrupt,  but 
the  ultimate  reduction  in  lift  and  lift-drag  ratio  is  comparable  to  that 
of  ventilated  flow.   Even  thin-airfoil  sections  of  small  design-lift 
coefficient  enter  this  cavltating  regime  of  poor  lift-drag  ratios  at 
speeds  in  excess  of  about  80  knots. 

Since  the  take-off  speed  of  modern  supersonic  aircraft  is  in  the 
range  of  I5O  to  200  knots,  lifting  surfaces  with  ca-vitating  or  venti- 
lating characteristics  superior  to  those  of  conventional  airfoil  sec- 
tions are  desirable.   Fortunately,  theoretical  work  by  Tulin  and  Burkart 
(ref .  1)  has  shown  that  superior  configurations  do  exist  and  they  have 
selected  a  cambered  configuration  for  operation  in  cavitating  or  venti- 
lated flow  which,  in  an  inviscid  fluid,  has  two-dimensional  lift-drag 
ratios  at  its  design  angle  of  attack  and  zero  cavitation  number,  about 
six  times  that  of  a  flat  plate.   If  such  a  cambered  hydrofoil  can  be 
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induced  to  ventilate  at  very  low  speeds,  while  the  aircraft  hull  still 
supports  most  of  the  load,  a  stable  and  efficient  take-off  run  may  be 
possible. 

The  present  report  is  concerned  with  an  experimental  and  theoretical 
investigation  of  one  of  the  Tiilin-Biirkart  cambered  configurations  and  a 
comparison  of  its  characteristics  with  those  of  a  flat  plate.   The  pur- 
pose of  the  investigation  was  (l)  to  find  out  if  the  superiority  of  the 
cambered  configurations  obtained  from  two-dimensional,  infinite -depth 
theory  would  be  maintained  at  finite  depth  and  aspect  ratio,  and  (2)  to 
develop  expressions  for  the  prediction  of  the  force  and  moment  character- 
istics of  a  ventilated  hydrofoil  of  arbitrary  section  profile  which 
woiild  include  these  three-dimensional  effects. 


SYMBOLS 


A  aspect  ratio 

An,AQ        coefficients  of  sine-series  expansion  of  vorticity  distri- 
bution on  equivalent  airfoil  section  at  infinite  depth, 
that  is, 

n(x)  =  2V(Aq  cot  ^  +  Ai  sin  9  +  A2  sin  29  ...  A^  sin  n9 


where 


A_  =  i  /   ^  cos  n9  d9 

^     rt   Jq    dx 


^     rt  Jn  dx 


0 

■'^i  h^-^  h     coefficients  of  sine-series  expansion  of  vorticity  distribu- 
tion on  hydrofoil  section 

a  distance  from  equivalent  airfoil  leading  edge  to  center  of 

pressure  in  chords 

Bj^,Bo        coefficients  of  cosine  series  defining  location  of  image 
vortex  in  airfoil  plane 

b  parameter  defining  location  of  spray  at  infinity  in  t  plane 

(see  ref.  2) 
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C-Q  total-drag  coefficient,  — 

Cj)  3_         profile-drag  coefficient  (total  drag  less  crossflow  and 
skin  friction),   -— 

Cj)  c         cross  flow -drag  coefficient  (normal  to  the  section) ,  — ^ 

Df 
Cf  skin-friction  drag  coefficient,  -^ 

Ct  total-lift  coefficient,  ■— 

Ct  total-lift  coefficient  of  equivalent  airfoil  section,  — 

Ct  1  lift  coefficient  exclusive  of  crossflow,  — = 

•L,  -L  q.S 

Ct  ^         crossflow  lift  coefficient,  -2. 
^y^  qS 

Cjn  pitching-moment  coefficient  (about  the  leading  edge),  

qSc 

Cm  pitching-moment  coefficient  of  equivalent  airfoil  section 

(about  the  leading  edge),  

qSc 

C^  ^  third-moment  coefficient  of  eqiiivalent  airfoil  section 


(see  ref.    l) ,     — 2_ 
qSc^ 

Cjj  resultant-force   coefficient  on  arbitrary  section,     -£- 

qS 

Cjj  f         resultant-force  coefficient  of  flat  plate,  -^ 

Cn,Co        coefficients  of  sine-series  expansion  of  vorticity  distri- 
bution on  equivalent  airfoil  section  at  arbitrary  depth 

c  chord 

D  total  drag  force 

Di  drag  force  excluding  crossflow  drag  and  skin  friction 
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Dq  drag  force  due  to  crossflow  (normal  to  section) 

Bf  drag  force  due  to  skin  friction 

d  leading-edge  depth  of  submersion 

E  Jones '  edge  correction,  ratio  of  seraiperimeter  to  span 
(see  ref.  Ik) 

F  resultant  force 

g  acceleration  due  to  gravity 

L  total  lift  force 

Lq  lift  force  due  to  crossflow 

L3_  lift  force  exclusive  of  crossflow,   L  -  L^ 

I  distance  from  hydrofoil  leading  edge  to  stagnation  point 
in  chords 

M  moment  about  leading  edge 

M2  third  moment  about  leading  edge,  2  /   p(x)x  dx 

m  =  CWa 

p  pressure,  Ib/sq  ft 

Pj,  pressure  within  cavity,  Ib/sq  ft 

Pq  pressure  at  mean  depth  of  hydrofoil,  Ib/sq  ft 

p  fluid  vapor  pressure,  Ib/sq  ft 

q  free-stream  dynamic  pressure,  -^pv 

r  distance  from  image  vortex  to  point  on  equivalent  airfoil 
in  chords 

S  area,  sq  ft 

s    _  span,  ft 

u  perturbation  velocity  in  X-dlrection 
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V  speed  of  advance,  fps 

V  perturbation  velocity  in  Y-direction 
X,Y  coordinate  axes 

X  distance  from  leading  edge  along  X-axis 

X  distance  from  leading  edge  to  center  of  pressiire 

y,  distance  from  the  X-axis  to  hydrofoil  lower  surface 

y^  distance  from  the  X-axis  to  hydrofoil  upper  surface 

a  geometric  angle  of  attack,  radians  unless  otherwise  specified 

ttg  angle-of -attack  increase  due  to  camber,      [Cq'   +  Cj^  -  C2/2J 

or  at  infinite  depth  0-^00=  (-^1  ~  ^/'^])   radians  unless 
otherwise  specified 

a^^  induced  angle  of  attack,  radians  unless  otherwise  specified 

cto  angle  between  hydrofoil  chord  line  and  reference  line,  posi- 

tive when  chord  line  is  below  reference  line,  rad-ians 
unless  otherwise  specified 

a'  angle  of  attack  measiired  from  hydrofoil  chord  line, 

a'  =  a  +  Oq,  radians  unless  otherwise  specified 

r  circulation,  strength  of  single  vortex 

7  central  angle  subtending  chord  of  circ\ilar-arc  hydrofoil 

6  spray  thickness  at  infinite  distance  downstream 

e  deviation  of  resultant-force  vector  from  normal  to  hydrofoil 

reference  line 


parameter  defining  distance  along  airfoil  chord, 
5  9) 

lb-sec^ 


X  =  -(1  -  cos  9) 
2 


mass  density. 


ft^ 

Po  -  Pc 


a  cavitation  number, 

q. 

a^  cavitation  n\miber  at   inception 
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T  correction  factor  for  variation  from  elliptical  plan  form 

(p  angle  between  spray  and  horizontal 

i|f  angle  between  x-axis  and  line  joining  image  vortex  with  a 

point  on  eqixLvalent  airfoil 

fl  vorticity 

I  >  indicates  "function  of"  for  example,   Cfj|a.[-  =  %  Aql   +  (xA 


Z  =  X  +  iY 

,  .   tan-l  £2 

=  a  +  e 

Subscripts : 

e 

effective 

0 

zero  depth  of  submersion 

t 

total 

CO 

infinite  depth  of  submersion 

c.p. 

center  of  pressijre 

Barred  symbols  refer  to  equivalent  airfoil  section  and  unbarred 
symbols  refer  to  the  supercavitating  hydrofoil  section. 


DESCRIPTION  OF  SUPERCAVITATING  FLOW 

Po  "  Pc 
The  parameter  defining  cavity  flow  is   ct  =  where  p^  is 

the  pressure  at  the  mean  depth,   p   the  pressure  within  the  cavity,  and 
q,  the  dynamic  pressure.   The  magnitude  of  a     for  the  condition  at  which 
cavitation  is  incipient  is  defined  by  the  particular  value  cr^.   If  c 
is  reduced  below  ff^,  cavitation  becomes  more  severe;  that  is,  the  cavita- 
tion zone  extends  over  a  larger  area.  When  a  hydrofoil  operates  at  suf- 
ficiently low  values  of  a,  the  cavity  formed  may  completely  enclose  the 
upper  or  suction  surface  and  extend  several  chords  downstream  as  shown  in 
figure  l(a) .   Theoretically  if  the  cavitation  parameter  is  reduced  to  zero, 
the  cavity  formed  will  extend  to  infinity.   The  flow  regime  where  the 
cavity  length  exceeds  the  chord  is  defined  as  supercavitating  flow. 
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Supercavitating  flow^  may  he   obtained  by  either  increasing  velocity 
or  cavity  pressure  or  both.   At  a  constant  depth  and  water  temperature, 
a  and  therefore,  the  length  of  the  cavity  is  dependent  only  on  the 
velocity  since  Pq  -  p^  is  then  Pq  -  p,^  and  is  constant.   If  part  or 
all  of  the  boundary  layer  of  a  configuration  is  separated,  the  eddying 
fluid  in  the  separated  region  can  be  replaced  by  a  continuo\is  flow  of 
lighter  fluid  such  as  air.   (See  refs.  5  and  k.)      Regulation  of  the  amount 
of  air  supplied  will  control  the  cavity  pressure  and  thus  the  length  of 
the  cavity  formed.   If  the  quantity  of  air  supplied  is  very  large,  the 
cavity  pressure  will  approach  the  ambient  pressure  Pq  and  a  very  long 
cavity  will  result  even  at  low-stream  velocities. 

The  ventilation  of  surface  piercing  hydrofoils  is  therefore  a  super- 
cavitating flow  due  to  large  quantities  of  air  supplied  from  the  atmos- 
phere to  separated  flow  on  the  suction  surface  of  the  hydrofoil.   Super- 
cavitating flow  as  a  result  of  ventilation  also  occurs  when  a  nonsxirface- 
piercing  hydrofoil  of  moderate  aspect  ratio  operates  near  the  free  surface 
(see  fig.  1(b)).   As  pointed  out  in  reference  5j  aii"  is  entrained  in  the 
trailing  vortices  and  drawn  to  the  suction  side  of  the  hydrofoil  causing 
a  long  trailing  cavity  to  enclose  completely  the  upper  surface  of  the 
hydrofoil  and  extend  far  downstream.   The  ventilated  type  cavity  described 
in  reference  5  differs  in  shape  from  those  formed  in  deeply  submerged  flow 
because  of  the  proximity  of  the  free  water  surface.   It  is  similar  to 
planing,  with  the  spray  forming  the  upper  surface  of  the  cavity.   Since 
the  cavity  pressure  is  approximately  the  same  as  the  ambient  pressure 
(at  small  depths  of  submersion),  the  cavitation  number  for  this  type 
flow  is  nearly  zero. 

The  experimental  portion  of  the  present  investigation  utilized  both 
of  the  methods  of  obtaining  low  values  of  a     discussed  in  the  foregoing 
paragraphs;  that  is,  varying  velocity  and  cavity  pressure.   However, 
the  case  of  zero  cavitation  number  was  of  particular  interest  since  it 
represents  the  condition  of  minimiJim  lift  and  also  because  a  major  portion 
of  the  existing  cavity  flow  theory  assumes  this  condition.   Therefore, 
most  of  this  report  treats  the  ventilated  type  supercavitating  flow. 


ZERO-CAVITATION-NUMBER  THEORY 


TWO-DIMENSIONAL  THEORY 


Flat  Plate 

The  characteristics  of  a  two-dimensional  inclined  flat  plate  in  an 
infinite  fluid,  operating  at  zero  cavitation  number,  have  been  obtained 
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by  Klrchoff  ajid  Rayleigh  (ref .  6) .   The  resultant  force  on  the  plate  is 
given  "by  the  well-known  equation 


I-,  2rt  sin  a 

CN,f  =  1 

'    4  +  Jt  sin  a 


(1) 


The  distance  from  the  leading  edge  to  the  center  of  pressure  is 


=  i(2 


c.p. 


3   cos 


U  +  jt  sin 


— ) 

n  a/ 


(2) 


Similar  work  was  performed  by  Green  (refs.  2  and  7)  to  include  the 
effect  of  the  free  water  surface  (but  neglecting  gravity) .   The  solution 
is  necessarily  obtained  in  terms  of  the  spray  thickness  6  rather  than 
the  more  useful  depth  of  submersion  and  is  given  as  two  parametric  equa- 
tions in  terms  of  the  parameter  b. 


Ct  =  C 


N,f 


cos  a 


2Vb 


(■ 


1/sin  a  cos  a 


=  ma 


(3a) 


S 
c 


where 

D 


.[^ 


b  -  l^sin  a  + 


b  -  cos  g 
D 


2  cos  a  +  (b  cos  a  -  l) logg 


(3t) 


b  +  1 


This  result  is  plotted  as  the  variation  of  m  with  5/c  for  various 
angles  of  attack  in  figure  2. 

Since  gravity  is  neglected  in  Green's  solution,  it  is  important  to 
understand  the  effect  of  the  absence  of  this  force  on  the  result  obtained. 
If  the  Froude  number  based  on  chord  V^/gc  is  large,  then  near  the  plate, 
where  the  streamline  radii  of  curvature  are  small,  the  inertia  forces  are 
indeed  large  compared  to  the  force  of  gravity.   Thus,  under  these  condi- 
tions, equations  (5)  for  the  force  on  the  plate  should  be  applicable,  for 
a  given  6/c,  in  spite  of  the  neglect  of  gravity  in  its  derivation.   How- 
ever, the  variable  6/c  is  not  usually  known.   It  is  desired  to  know  the 
force  on  the  plate  for  a  given  depth  of  submersion;  therefore,  the  rela- 
tionship between  the  leading-edge  depth  of  submersion  and  the  spray 
thickness  is  needed.   The  effect  of  gravity  becomes  very  important  in 
determining  this  relationship.   Because  of  the  neglect  of  gravity  in  the 
analysis  of  Green,  the  plate  is  always  located  infinitely  above  the  still 
water  surface  -  a  physical  impossibility  in  a  real  fluid.   Gravity  cannot 
be  neglected  in  determining  the  location  of  the  free  water  surface.   The 
reason  for  this  is  that  several  chords  away  from  the  plate  the  streamline 
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raxill  of  ciirvatixre  become  very  large  and  thus  the  inertia  forces  become 
of  the  same  or  lower  order  of  magnitude  than  the  force  of  gravity.   When 
gravity  is  present  the  distance  of  the  plate  above  the  water  surface  at 
infinity  is  at  least  limited  to  the  dynamic  pressure  head.   The  effect 
of  finite  aspect  ratio  is  also  to  place  the  plate  finitely  near  the  free 
water  surface  as  will  be  discussed  in  a  subsequent  section. 

In  summary,  if  the  Froude  number  is  large,  the  forces  on  an  inclined 
flat  plate  near  a  free  water  surface  can  be  obtained  in  terms  of  the  spray 
thickness   6/c  from  equations  (3) .   The  relationship  between  the  spray 
thickness  and  the  actual  leading-edge  depth  of  submersion  cannot  be  deter- 
mined from  Green's  analysis,  but  the  relationship  is  known  to  be  influ- 
enced by  the  presence  of  gravity  and  the  aspect  ratio  of  the  plate.   The 
method  by  which  the  lines  of  constant  depth  (shown  in  fig.  2)  were  obtained 
is  discussed  in  a  subsequent  section. 

The  angle  0  that  the  spray  makes  with  the  direction  of  motion  is 
also  obtained  from  reference  2  as 


d   =  cos-^f^  ^°^  '^  -  A  -   a 
\  b  -  cos  a  / 


(M 


The  variation  of  <ji      (eq.  (1+))  with  6/c   (eq.  (jb))  is  shown  in  figure  3. 


Cambered  Sections 

Cambered  surfaces  can  theoretically  be  analyzed  in  two  dimensions  by 
the  method  of  Levi-Civita  (ref .  7) •   However,  like  many  conformal  mapping 
problems  the  method  is  very  difficult  to  apply  to  a  particular  configura- 
tion and  only  a  few  specific  solutions  have  been  obtained.   Among  these 
is  the  work  of  Rosenhead  (ref.  8)  and  Wu  (ref.  9)-   Although  the  solution 
of  Wu  is  applicable  in  principle  to  arbitrary  sections,  the  solution  was 
presented  only  for  the  circular  arc.   A  particular  advantage  of  Wu's 
solution  is  that  it  includes  the  effects  of  nonzero  cavitation  number. 

The  most  useful  treatment  of  cambered  s\arfaces  is  the  linearized 
theory  of  Tulin  and  Burkart  (ref.  1)  which  is  readily  applicable  to  any 
surface  configuration  (with  positive  lower  surface  pressures)  as  long 
as  the  angle  of  attack  and  camber  are  small.   The  principal  results  of 
this  linearized  theory  are  summarized  below. 

The  supercavitating  hydrofoil  problem  in  the  Z-plane  is_transformed 
into  an  airfoil  problem  in  the  Z-plajie  by  the  relationship  Z  =  -/Z.   By 
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denoting  properties  of  the  equivalent  airfoil  with  barred  symbols  and 
those  of  the  hydrofoil  with  unbarred  symbols,  the  following  relationships 
are  derived: 


u(x)  =  u(x2)  (6) 

Cl=  c„j  =  |(ao  +  Ai -:^)  (7) 


^^-U'^f-i{^*¥f 


(8) 


Cm  =  Cni,3  =  32?^   ^  ^^^  "  '^^^   "^  ^^5  "  2)  ^^^ 


The  coefficients  A^^  are  the  thin-airfoil  coefficients  in  the  sine- 
series  expansion  of  the  airfoil-vorticity  distribution 


Q(x)   =  2V(Aq  cot  I  +  \  A^   sin  n9J  (lO) 


n=l 

where 


X  =  |(1  -  cos   9)  (0  ^  9  ^  rt) 

and  can  be  foimd  for  a  given  configuration  from  the  following  eq\iations 

AQ  =  -ir      ^d9+a  =  a  +  AQ'  (lla) 


12 
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An  -  rt 


-^  cos  n9  de 
dx 


(lib) 


By  -using  equations  (7)  and  (8),  the  condition  that  Aq  in  equation  (lO) 
is  zero,  that  Q,{x)      does  not  contain  harmonics  greater  than  two,  and 
that  Q,{x)      is  everywhere  positive  over  the  chord  (that  is,  positive 
lower  s\irface  pressures  on  the  hydrofoil),  Tulin  and  Burkart  (ref.  l) 
obtained  a  low-drag  family  of  hydrofoils.   This  particular  family  of 
hydrofoils  is  given  by  the  equation 


y  _  Ai 


"lie 


x\5/2 


Mf 


(12) 


The  lift-drag  ratio  of  this  foil  at  its  design  incidence  (that  is,  a  =  O) 
is  from  equations  (7)  and  (8) 


L  ^  2-?/  fl 


h    2Ct 


(13) 


This  6.25-fold  improvement  over  the  L/D  of  a  flat  plate  is  most  encour- 
aging and  was  the  impetus  for  this  investigation. 

The  restrictions  imposed  by  the  assumptions  of  the  T\ilin-Burkart 
theory  prevents  its  use  in  the  calculation  of  the  characteristics  of 
hydrofoils  siiitable  for  use  as  aircraft  landing  gear.   Here,  because  of 
the  high  hydrofoil  loads  on  necessarily  thin  hydrofoils  the  aspect  ratio 
may  be  as  low  as  1  or  2.   Also  the  hydrofoil  must  operate  near  the  free 
water  surface  and  in  some  instances  at  large  angles  of  attack.   Thus, 
the  effects  of  these  variables  on  the  characteristics  of  supercavitating 
hydrofoils  (particularly  of  cambered  sections)  is  needed.   Much  of  this 
information  can  be  obtained  by  certain  modifications  to  the  existing  two- 
dimensional  theory  discussed  in  preceding  paragraphs. 


MODIFICATIONS  OF  TWO-DIMENSIONAL  THEORY 


Nonlinear  Equation  for  Lift  at  Infinite  Depth 


For  any  configuration  the  reference  line  of  the  section  from  which 
the  angle  of  attack  is  measured  can  be  chosen  such  that  Aq'  =  0.   It  is 
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convenient  in  this  report  to  assume  that  a  =  0  refers  to  the  reference 
line  which  makes  Aq'  =  0.  With  this  assumption,  equation  (7)  may  be 
written  as 


cl  =  if  ^  Ai  -  ^)  =  f  (^  +  -c 


(Ik) 


where  a^  is  the  effective  increase  in  angle  of  attack  due  to  camber 


(^^  -  *) 


Thus,  the  solution  for  cambered  hydrofoils  is  merely  the  flat 

plate  linearized  solution  -^a  with  a  replaced  by  a  +  a_.   This  is 

exactly  analogous  to  the  influence  of  camber  on  airfoils  in  an  infinite 
fluid  where  there  is  an  effective  increase  in  angle  of  attack  due  to  the 
camber.   By  carrying  this  procedure  further,  and  by  applying  it  to  the 
resioltant  force  rather  than  the  lift,  the  nonlinear  solution  of  Ray lei gh 
becomes  applicable  to  arbitrary  configurations  simply  by  replacing  a 
by  a  +  a^,;  that  is, 

2n  sin(a  +  a^) 

^N  -  I 7^7 -\  ^15) 

4  +  jt  sin(a  +  a,,) 


The  lift  will  then  be 


^   2Jt  sin(a  +  oc)  ,  ,. 

Cl  =  cos  3  (16) 

k  +  It   sin(a  +  ac) 


Here,   p  =  a  +  e,  where  e  denotes  the  deviation  of  the  resultant-force 
vector  from  the  normal  to  the  hydrofoil  reference  line.   For  large  values 
of  a,   e  is  small  compared  with  a  and  cos  p  =  cos  a.  When  a  is 
very  small,  about  0°,   €  is  a  maxiraijm  and  can  be  shown  by  use  of  equa- 

A-|2 
tions  (7)  and  (8)  to  be  — - — = r-  if  Aq'  =  0  as  assiomed.   For  any 

K*x  -  ¥) 

practical  low-drag  hydrofoil,  this  value  will  almost  always  be  less  than 
5°  for  which  the  cosine  is  very  nearly  1  or  cosine  (a  +  e)  =  cos  a  =  1. 
Therefore,   cos  p  in  equation  (I6)  may  be  replaced  by  cos  a  with  little 
loss  in  accuracy  and  great  gain  in  simplicity.   Equation  (16)  then  becomes 

^   2it  sinfa  +  clq)  /   \ 

Cl  -  : -, ^  cos  a  (17) 

4  +  Jt  sin(a  +  ac) 
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For  a  circular-arc  hydrofoil  of  central  angle  7,  it  is  shown  in 
appendix  A  that  ac  =  (9/16)7  and  that  the  reference  line  must  be  chosen 
at  an  angle  7/8  to  the  chord  line  so  that  Aq'  =  0.   The  result  obtained 
by  substituting  this  value  into  equation  (I7)  is  compared  in  figure  h 
with  the  linear  solution  of  Tulin  and  Burkart  (eq.  (7))  and  the  nonlinear 
solution  of  Wu  (ref.  9)  for  two  circular-arc  profiles.  The  agreement  of 
equation  (I7)  with  the  more  exact  solution  of  Wu  is  good  over  the  entire 
range  of  angle  of  attack  from  0°  to  90° •   Similar  agreement  is  expected 
for  any  configuration  of  small  camber. 


Nonlinear  Equation  For  Lift  at  Finite  Depth 

The  successful  modification  of  the  Rayleigh  equation  to  include 
cambered  configurations  leads  at  once  to  a  similar  modification  of  the 
solution  of  Green.   However,  in  this  case  the  argument  for  replacing  a 
by  (a  +  a^)   is  very  weak  unless  the  section  coefficients  which  deter- 
mine ac  are  known  as  a  function  of  the  depth  of  submersion. 

An  examination  of  the  linearized  expressions  for  the  lift  coefficient 
of  arbitrary  foils  at  infinite  depth  and  at  zero  depth  reveals  that  both 
the  lift-curve  slope  and  the  increase  in  angle  of  attack  due  to  camber 
do  change  with  depth  of  submersion.   At  infinite  depth  the  linearized 
expression  for  lift  coefficient  is  given  by  equation  (7) •   At  zero  depth 
the  lift  coefficient  must  be  one  half  of  the  fully  wetted  value  obtained 
from  thin-airfoil  theory  as  pointed  out  in  reference  10;  that  is. 


Cl,0  =  «  Ao,h  +^  (18) 


where  Aq  5^  and  A^  ^  ^^^  "the  thin-airfoil  coefficients  of  the  section 
in  the  hydrofoil  plane  and  are  given  by  the  expressions 


Ao,h 


-n   X'  i  ^^  (^9a) 


For  the  T;illn- Burkart  section  at  zero  ajigle  of  attack  these  values  may 
be  determined  as 
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Aq^^  =  0.227A^  (20a) 

Ai^h  =  I.I5IA-L  (20b) 

Thus,  from  equations  (7)  and  (I8)  it  is  seen  that,  for  a  flat  plate  at 
small  angles,  the  lift  coefficient  goes  from  ^  at  infinite  depth  to 
ncL     at  zero  depth  (as  given  by  Green)  whereas,  for  the  Tulin-Burkart 
section  at  zero  angle  of  attack  these  values  are  •2(l.25Ai)   at  infinite 
depth  and  «(0.802Ai)   at  zero  depth.   Although  the  flat-plate  lift  coef- 
ficient doubles  in  going  from  infinite  to  zero  depth,  the  ratio  is  only 
1.28  for  the  cambered  section.   The  important  point  to  note  is  that  the 
value  of  ac  for  the  Tulin-Burkart  section  changes  from  1.25A2_  to 
0.802A-L. 

It  is  now  desirable  to  determine  a^,  for  finite  depths  of  submersion. 
This  can  be  accomplished  by  modifying  the  linearized  theory  of  reference  1 
to  include  the  effects  of  the  free  water  surface. 

The  effect  of  the  free  water  surface  may  be  obtained  by  finding  the 
transformation  which_will  map  the  free  water  surface,  the  hydrofoil,  and 
the  cavity  into  the  X-axis  of  the  Z-plane.   The  required  transformation 
is  of  the  form  Z  =  Z  +  K-]_  logg  Z  +  K^.   By  using  this  transformation  it 
is  possible,  in  principle,  to  determine  the  forces  on  the  hydrofoil  by 
following  the  procedure  used  in  reference  1.   However,  the  solution  by 
this  method  became  very  cumbersome  and  was  abandoned  when  an  approximate 
method  was  discovered.   The  approximate  method  continues  with  the  trans- 
formation Z  =  -|z  used  in  reference  1. 

In  figure  5  it  "lay  be  seen  that  Z  =  -/z  transforms  the  free  water 
surface  in  the  hydrofoil  plane,  where  u  =  0  (see  fig.  5(a)),  into  a 
hyperbola  in  the  third  quadrant  of  the  airfoil  plane  (fig.  '^{'b))  .      The 
boundary  condition  that  must  be  satisfied  on  this  hyperbola  is  that  the 
perturbation  velocity  u  be  zero  because  in  complex  velocity  problems 
of  the  type  considered  here,  it  is  the  lines  of  constant  velocity  which 
are  being  transformed  and  not  lines  of  constaxit  velocity  potential  or 
stream  fiinction. 

For  the  particular  case  of  zero  depth  the  hydrofoil  problem  is  trans- 
formed by  Z  =  -Jz  into  the  fourth  quadrant  of  the  airfoil  plane.   Thus, 
it  may  be  seen  in  figure  5(b)  that  the  free  water  surface  adds  the  con- 
dition that  u  =  0  along  the  negative  Y-a:xis.   This  additional  boundary 
condition  csui  be  satisfied  (along  with  the  other  infinite-depth  boundary 
conditions)  by  locating  an  image  of  the  airf oil- vorti city  distribution 
along  the  negative  X-ajcis.   The  direction  of  this  vorticity  must  be 
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opposite  to  that  of  the  airfoil  in  order  to  make  u  =  0  at  all  points 
along  the  negative  Y-axis.   A  simpler  and  often  used  approximation  is  to 
replace  the  distributed  image  vorticity  by  a  single  vortex,  eqiial  in 
strength  to  the  airfoil  circiilation,  at  a  location  equal  to  the  distance 
from  the  leading  edge  to  the  airfoil  center  of  pressure  as  shown  in 
figure  5(b) . 

For  finite  depth  of  submersion,  the  condition  that  u  =  0  must  be 
satisfied  at  all  points  on  the  hyperbola  and  on  the  negative  X-axis.   It 
is  not  possible  to  satisfy  these  conditions  with  a  single  vortex  as  was 
done  for  the  case  of  zero  depth.   However,  the  influence  on  the  airfoil 
of  the  infinite  array  of  vortices  needed  to  satisfy  the  boundary  condi- 
tions shown  in  figiore  6  may  be  approximated  by  a  single  vortex  of  strength 
r  in  the  location  shown.   The  adequacy  of  the  approximation  can  be  deter- 
mined by  calculating  the  effect  of  this  image  vortex  on  the  lift  of  a 
flat  plate  as  the  depth  of  submersion  is  varied  and  then,  comparing  the 
result  with  the  exact  solution  of  Green.   This  depth  effect  may  be  deter- 
mined by  concentrating  the  airfoil  circulation  at  its  center  of  pressure 

and  the  image  circulation  at  a  point  ^  +  cj—  forward  of  the  leading 

~  fd  4    »c 

edge  and  •^j—  below  the  leading  edge  (see  fig.  6)  and  computing  the 

''  ■'in 

total  downwash  on  the  flat  plate  at  its  ^—   point  (see  ref .  11) .   The 

4 

method  assumes  that  the  center  of  pressure  of  the  airfoil  remains  con- 
stant at  c/k     as  the  depth  changes.  The  resulting  downwash  angle  a 

at  the  ^  point  is  calculated  to  be 


itcV 


1-i/ 


1  +  /d/c 


+  2^  +  (5A)(d/c) 


(21) 


Since  T^  =   ncaV,  the  ratio  of  hydrofoil  lift  at  finite  depth  Cl  to  the 
lift  at  infinite  depth  is 


•"-m  V  L 


T._oo  ^m  _  00  -^n-^  ^  QQ 


^L,oo  Om,oo  ^^^^ 


(22) 


therefore 


Cl 


Cl,oo      r„      ^  ^  ^/ 1  +  y/STc ) 

2\i  +  2/d7^  +  (5A)(d/c)y 


(23) 
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Equation  (23)  is  compared  in  figijre  7  with  the  exact  solution  of  Green 
(see  fig.  2.,     a  =  O)  ,  and  the  agreement  is  excellent. 

The  adequacy  of  the  method  used  in  determining  the  influence  of  free 
water  surface  proximity  on  the  lift  of  a  flat  plate  justifies  its  use  on 
cambered  foils.   However,  for  cambered  foils  the  problem  is  more  difficult 

because  a  control  point  such  as  '^—  for  a  flat  plate  is  not  sufficient  to 

4 
determine  the  section  coefficients  defining  the  vorticity  distribution. 

This  final  vorticity  distribution  n(x)   resulting  from  the  camber  and  the 
image  vortex  must  be  determined,  particularly  if  a  knowledge  of  the  pres- 
sure distribution  and  thus  the  drag  and  center  of  pressure  is  desired. 

In  figure  8,  it  may  be  assumed  that  the  final  vorticity  distribution 
on  the  equivalent  airfoil  is  given  by  the  equation 

fl(x)  =  2V  Cq  cot  i  +  \  Q,^   sin  n9|  {2h) 

Then  the  induced  velocity  at  x'   due  to  fi(x)   is,  from  thin-airfoil 
theory  (see  ref.  12), 

V  =  vl-Cn  +  ^   C^  cos  nel  (25) 


^0  -  ^  -n 
n=l 


The  resulting  total  circulation  due  to  the  vorticity  given  by  equation  (2^+) 
is,  from  thin-airfoil  theory, 

(26) 


Since  the  image  vortex  has  a  strength  equal  and  opposite  to  P,  the 
velocity  induced  by  rj[  at  a  point  x'   on  the  airfoil  is 

v-t  =  cos  ilr  (27) 

^       2jrrc 
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vhere 


or 


where 


or 


=  i/M!^if^) 


a  4.  M  +  41 

I  c    c 

COS  \(f  =  


COS  i{f         V  c 


a  +.4  +4^ 


„  ,  ,d  ^  X'   ^  1  d 


(28) 


Cp 
Co  +  Ci  -  ^ 

a  =  ^ ^  (29) 

Co  +^ 


Replacing  x'      by  S.(i  _  cos  9),  equation  (28)  can  be  expanded  in  a 
Fourier  series  as 

CO 

H£|J  =  f(9)  =  Bo  +  ^  Bn  cos  ne  (50) 

n=l 


where 


Bo  =  i  JJ  f(9)d9 
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and 


B  =  -  /   f(e)cos  n9  de 

^  Jo 


By  sulDstltuting  equations  (26)  and  (jO)  into  equation  (27)  gives 


VICq  + 


VI  =  — i-- ^(Bo  +  )  Bn  cos  n9l  (5I) 


n=l 


The  equivalent  airfoil  slope  -^  when  expanded  in  a  cosine  series 
is  ^ 

00 
If  =  -Aq  +  y  An  c°s  ne  (52) 


dx         _ 
n=l 


Equating  the  resulting  streamline  slope  to  the  equivalent  airfoil  slope 
gives 


It  ^  UJl  ^  d^  (33) 

V     V     dx 


The  substitution  of  equations  (25),  (51)^  an<i  (52)  into  equation  (53) 
gives 


-Co  +  )  Cn  cos  ne  +  l \Bq   +  >  Bn  cos  n9  =  -Aq   +  >  An  cos  nS 

n=l  \    n=l        /        n=l 

By  equating  coefficients  of  like  terms,  the  C  coefficients  are  deter- 
mined as 

^     1^  +  B^  -  2Bq 
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i^  +  B-L  -  2Bo 


'-  =  ^n  -  ^^   '1  "  (35c) 

U  +  B-L  -  2Bo 

If  Aq  =  a  and  A^  =  0  where  n  ^  1,  the  effect  of  depth  of  sub- 
mersion on  the  lift  coefficient  of  a  flat  plate  can  be  computed  from  the 
coefficients  obtained  from  equations  (55).   The  values  of        computed 

by  this  method  have  been  found  to  be  in  excellent  agreement  with  the  solu- 
tion of  Green  for  a  =  0. 

For  the  particular  condition  of  Aq  =  0   (the  case  of  hydrofoils  such 
as  the  Tulin-Burkart  section  at  zero  angle  of  attack)  equations  (55) 
become 

BdAt 

Co  = ^-^ (56a) 

I4-  +  B-L  -  2Bq 


2A-,(2  -  Bq)  ,  ,  , 


i|  +  Bn  -  2Br 


Cn  =  An  - ^^ (56c) 


ij-  +  BjL  -  2B(; 


The  coefficients  Bq  and  Bj^  as  obtained  from  equation  (50)  are  plotted 
in  figure  9  against  the  distance  to  the  center  of  pressiire  a  for  several 
depths  of  submersion.   For  the  special  condition  of  zero  depth  the  B 
coefficients  for  the  Tulln-Burkart  section  are  found  to  be  I.296,  O.772, 
and  0.22  for  Bq,   B-]_,  and  B^,   respectively.   In  making  this  computation, 
the  final  center-of-pressirre  location  a  is  used;  therefore,   a  is  first 
given  an  assumed  value,  the  B  and  C  coefficients  determined,  then 
from  the  resulting  C  coefficients,   a  is  calculated  from  equation  (29) 
and  the  procedure  repeated  if  necessary.  For  the  Tulin-Burkaxt  section 
a  is  found  to  be  0.U2  for  d/c  =0.   By  using  the  final  B  coefficients 
the  C  coefficients  are  determined  from  equations  (56)  as  0.595A2_, 
0.6i4-6A-|^,  and  -O.6IA1  for  Cq,   C^,  and  C2,  respectively. 
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From  equation  (7)  the  ratio  of  the  Tulin-Burkart  section  lift  at 
zero  depth  to  the  lift  at  infinite  depth  is 

Cl,o   (0.595  +  0.646  +  0.305) Ai   ^  ,         ,  , 

C  ~- ^^ '  "■"'  *''' 


The  value  1.24  compares  favorably  with  the  more  exact  value  of  1.28  given 
in  the  paragraph  following  equations  (20) .   The  results  of  calculating 

for  finite  depths  of  submersion  for  the  Tulin-Burkart  section  are 

CL,oo 

plotted  in  figure  10.   Similar  calculations  have  been  made  for  a  circular 
arc,  and  the  configurations  given  in  reference  I3.   The  results  of  these 
calculations  are  also  shown  in  figure  10. 

The  true  linearized  lift-curve  slope  m  for  finite  depths  of  sub- 
mersion in  the  equation  C^  =  m(a-  +  a^)   is  that  shown  in  figure  2  for 
a  =  0.   Therefore,  the  effective  angle  of  attack  due  to  camber  a^,  is 
obtained  from  the  following  relationship: 

^c^O)°^c  _  ^L 
^  a,    '  %o. 

2   C,oo        ' 


Therefore 


"-c       rt     Cl 


o-c^oo   2m(a^o)  ^l^oo 


(38) 


(X, 


Values  of  are  plotted  against  d/c  in  figure  11  for  the  Tulin- 

Burkart,  the  circular  arc,  and  the  sections  given  in  reference  I3. 

Equation  (38)  is  obviously  limited  by  the  linearizing  assumptions 
made  in  its  derivation.   An  important  limitation  is  due  to  the  assumption 
that  the  free  surface  is  always  horizontal  and  thus   6/c  =  d/c.   At 
small  depth-chord  ratios  and  particularly  for  large  magnitudes  of  camber 
the  free  water  surface  is  not  horizontal  and  6/c  /  d/c.   Thus,  for  small 
values  of  d/c  and  large  magnitudes  of  camber  the  values  of  clq/o^q   m 
given  in  figure  11  are  probably  too  low. 

With  a  knowledge  of  the  angle  of  attack  due  to  camber  a^,  at  finite 
depths  of  submersion.  Green's  solution  is  now  modified  to  include  camber 
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by  treating  the  effective  angle  of  attack  as   (a  +  a^) ,  where  a^  is 
ohtained  from  figure  11.   This  is  exactly  the  method  used  in  modifying 
the  Rayleigh  equation  to  obtain  the  nonlinear  approximation  for  the  lift 
coefficient  at  infinite  depth.   With  this  assumption,  the  resultant-force 
coefficient  for  a  cambered  hydrofoil  at  any  positive  depth  of  submersion 
is  obtained  in  terms  of  the  spray  thickness   6/c  from  equations  (5)  as 


Cjjja}  ^  Cu^fja  +  acj  (39) 


Equation  (39)  states  that  the  resultant  force  on  a  cambered  section  is 
approximated  by  replacing  a  in  Green's  solution  for  a  flat  plate  by 
the  effective  angle  of  attack  (a  +  a^) •   It  will  be  shown  that  the 
resultant  force  will  deviate  only  slightly  from  the  normal  (as  previously 
pointed  out  for  the  condition  of  infinite  depth)  and  therefore 


Actually,  by  itself,  equation  (i+O)  is  of  little  practical  value 
because  C|^  f     is  given  in  equations  (3)  in  terms  of  the  usually  \inknawn 
parameter  5/c  instead  of  the  more  useful  d/c.   Thus,  it  is  necessary 
to  determine  the  relationship  between  6/c  and  d/c   or  at  least  an 
empirical  substitute. 


Relation  Between  Depth  of  Submersion  and  Spray  Thickness 

Certain  relationships  between  the  leading-edge  depth  of  submersion 
and  the  spray  thickness  of  a  flat  plate  can  be  stated  and  are  given 
as  follows : 

(1)  The  trivial  but  useful  case  for  a  =  0  where  d/c  =  6/c 

(2)  The  case  of  a  =  90°,   6/c  =  00,  where  the  stagnation  streamline 
is  parallel  to  the  direction  of  motion  and  in  order  to  satisfy  the  con- 
tinuity equation  6/c  =  d/c  +  O.5   (see  fig.  12(a)) 

(3)  The  case  of  a  =  90°,   6/c  =  0  where  again  the  stagnation  line 
(this  time  the  free  water  surface)  is  parallel  to  the  direction  of  motion 
and  &/c  =  d/c  +  1  (see  fig.  12(b)) 

Cases  (2)  and  (3)  suggest  that  the  stagnation  line  for  a  =  90°  is 
also  parallel  to  the  direction  of  motion  for  finite   6/c   (as  shown  by 
the  solid  line  in  fig.  12(c)).   However,  this  is  not  consistent  with  the 
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results  of  Green's  analysis  where  it  can  be  shovn  from  momentum  considera- 
tions that  the  streamline  is  curved  as  shovm  "by  the  broken  line  in  fig- 
ure 12(c).   However,  when  the  effects  of  gravity  and  even  slight  devia- 
tions from  infinite  aspect  ratio  are  considered,  the  true  location  of  the 
stagnation  line  must  lie  between  these  two  lines.   The  tendency  of  gravity 
to  cause  this  change  has  been  discussed  already.   Possibly  more  important 
than  gravity  is  the  influence  of  finite  aspect  ratio. 

If  the  supercavitating  hydrofoil  is  replaced  by  a  system  of  horseshoe 
vortices,  it  can  be  shown  that  (for  finite  angles  of  attack  and  infinite 
aspect  ratio)  the  stagnation  line  infinitely  forward  of  the  hydrofoil  is 
also  infinitely  below  the  stagnation  point  -  as  obtained  by  Green.   How- 
ever, for  an  aspect  ratio  of  as  much  as  100,  the  stagnation  line  approaches 
a  finite  asymptote  about  0.2  chord  below  the  stagnation  point  when  the 
lift  coefficient  is  about  O.5.   At  a  lift  coefficient  of  O.5,  the  asymp- 
tote becomes  about  0.1  chord  for  moderate  aspect  ratios.   This  was  observed 
in  the  experiments  described  in  a  subsequent  section  where  it  was  noted 
for  the  angle-of-attack  range  investigated  that  the  free  water  surface 
was  practically  undisturbed  only  slightly  forward  of  the  hydrofoil  leading 
edge.   Therefore  for  an  angle  of  attack  of  90°  where  the  total  circulation 
is  zero,  it  is  not  unreasonable  to  expect  the  solid  stagnation  line  in 
figure  12(c)  to  be  more  nearly  correct  than  the  dashed  one  when  both  finite 
aspect  ratio  and  gravity  are  present.   It  should  be  noted  that  although 
the  stagnation  line  may  de-viate  from  the  theory  as  given  by  Green,  that 
the  pressure  distribution  and  force  on  the  plate  are  assumed  to  remain 
unaltered  from  the  expressions  given  in  equations  (5) •   Therefore,  if  at 
90°  the  stagnation  line  is  assiomed  parallel  to  the  direction  of  motion. 
It  follows  that  for  a  =  90° 

6/c  =  d/c  +  Z/c  (i+1) 

where   z/c  is  the  dimensionless  distance  from  the  leading  edge  to  the' 
stagnation  line  (see  fig.  12(c)).   This  distance  can  be  obtained  from 
Green's  work  and  is,  for  a,  =  90° 

z/c  =  ^ \:^^ V^^ W         (^2) 

.(b  -  WT-^)  .  iog,(^) 

where  b  has  been  previously  defined  in  equations  (3) . 


1  +^ 
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Case  (1)  and  equations  (hi)    and  {k2)    give  the  end  points  (at  a  =  0° 
and  90°)  of  the  lines  of  constant  d/c  which  are  superimposed  on  fig- 
ure  2.      The  trends  of  these  cwrves   at  intermediate  angles  of  attack  were 
obtained  from  experiment  and  will  be  discussed  in  a  subsequent  section. 

The  preceding  discussion  in  this  report  has  been  concerned  with  the 
effects  of  gravity  and  finite  aspect  ratio  on  the  relationship  between  the 
spray  thickness  and  the  leading-edge  depth  of  submersion.   It  is  also  of 
interest  to  examine  at  least  qualitatively  the  effects  of  gravity  and 
finite  aspect  ratio  on  the  spray  angle  given  by  equation  (k)    and  plotted 
in  figure  3.   In  Green's  analysis,  it  is  shown  that  at  a  =  90°  the  total 
force  on  the  entire  plate  is  F  =  p5v  (l  -  cos  ^) ;  that  is,  the  force  is 
exactly  eqiial  to  the  change  in  momentum  due  only  to  the  flow  in  the  spray. 
This  means  that  the  stagnation  line  must  be  curved  as  pointed  out  in  the 
preceding  discussion.   The  straight  stagnation  line  (which  leads  to  equa- 
tion (Ul))  Implies  that  the  change  in  momentum  in  the  spray  flow  results 
only  in  the  force  on  that  portion  of  the  plate  above  the  stagnation 
point.   Since  the  change  in  momentum  of  the  spray  flow  now  represents 
only  a  portion  of  the  total  force  on  the  plate,  and  since  5/c  is  con- 
sidered to  be  the  same  for  both  stagnation-line  curvatures,  the  angle  ^ 
must  decrease  from  the  value  given  by  Green  in  figure  5.   That  is,  the 
angle  0  will  change  so  that  only  the  force  on  the  plate  above  the  stag- 
nation line  is  given  by  the  expression  F  =  p6V^(l  -  cos  0). 


THREE  DIMENSIONAL  THEORY  AT  FINITE  DEPTH 

Lift 

The  flow  about  a  supercavitating  hydrofoil  may  be  constructed  by  a 
suitable  combination  of  sources  and  vortices.  The  vortices  contribute 
unsymmetrical  velocity  components  and  lift,  the  sources  contribute 
symmetrical  components  which  provide  thickness  for  the  cavity  but  no 
lift.  For  a  finite  span  the  vortices  can  not  end  at  the  tips  of  the 
foil  and  a  system  of  horseshoe  vortices  must  be  combined  with  the  sources 
to  describe  the  flow.  If  it  is  assumed  that  the  Influence  of  finite  span 
on  the  two  dimensional  lift  coefficient  is  due  to  the  effects  of  the 
trailing  vortlcity  then  the  resulting  effect  of  aspect  ratio  is  exactly 
the  same  as  for  a  fully  wetted  airfoil.  Jones  (ref .  Ik)   gives  the  lift 
of  a  fully  wetted  elliptical  flat  plate  as 

C L  =  I  2n(a  -  a^)  (45) 

where  E  is  the  ratio  of  semiperimeter  to  the  span  and  aj_     is  the 

induced  angle  of  attack  caused  by  the  trailing  vortlcity.   Thus  the  effect 
of  aspect  ratio  is  to  decrease  the  two  dimensional  lift  curve  slope  by  a 
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factor  i/e  and  to  decrease  the  effective  angle  of  attack  by  an  incre- 
ment aj_.  Therefore  for  the  finite  aspect  ratio  supercavltatlng  hydrofoil 
at  infinite  depth,  equation  (ik)   is  modified  to  give 

or  more  generally  for  finite  depth,  equation  (kO)   becomes 

^L,l  =  I  ^Mi°'   +  ^c  -  a-i)cos  a  (45) 

where  for  rectangular  plan  form  of  aspect  ratio  A,  E  =  ^  "*"  -^  and 

a^   =  ^(1  +  t)  (46) 

^    itA 

where  t  is  a  correction  for  plan  form  (see  ref.  12). 

Another  effect  due  to  finite  aspect  ratio  is  the  concept  of  additional 
lift  due  to  crossflow  (see  refs.  I5  and  16).  This  crossflow  lift  is 
assumed  due  to  the  drag  on  the  hydrofoil  contributed  by  the  component  of 
free -stream  velocity  normal  to  the  hydrofoil.  In  the  present  case  of  zero 
cavitation  number,  the  crossflow  drag  coefficient  is  the  Rayleigh  value, 
0.88.  Since  this  lift  is  caused  only  by  the  spanwise  flow  (flow  around 
the  ends  of  the  plate)  it  is  also  modified  to  account  for  the  aspect  ratio 
by  the  Jones'  edge  correction,   i/e.  Since  only  the  spanwise  flow  is  con- 
sidered, E  is  now  the  ratio  of  semiperlmeter  to  chord.   Because  the  flow 
being  considered  is  normal  to  the  plate,  the  induced  angle  for  this  flow  is 
zero.  Thus  for  a  flat  plate,  the  crossflow  lift  C^  q     is 


C 


L,c 


i —  0,88  sln^  a  cos  a  (47) 


A  +  1 


No  experimental  or  theoretical  information  on  the  crossflow  lift  of 
cambered  surfaces  is  available  in  the  literature.   In  order  to  approximate 
this  component  the  following  assumptions  are  made: 

(1)  The  crossflow  force  acts  normal  to  the  hydrofoil  chord  line. 

(2)  The  effective  direction  of  the  free  stream  on  the  plate  is  altered 
by  the  increase  in  angle  of  attack  due  to  camber  a^ . 

Thus,  the  crossflow  lift  on  cambered  sections  is  assumed  to  be 
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Ct   „  =       -*-       0.88  sin2(a'    +  ac'Jcos  a'  (1+8) 

L,c       A  +  1  ^  ^J 

where  a'  =  a  +  aQ,  and  aQ  is  the  inclination  of  the  chord  line  to  the 
reference  line  of  the  section  (positive  if  the  chord  line  is  below  the 
reference  line)  and  a^  is  obtained  from  figure  11  for  the  depth  of 
interest. 

The  total  lift  on  a  finite  aspect  ratio  hydrofoil  operating  near  the 
free  water  surface  is  therefore  obtained  by  adding  equation  (1+8)  to  equa- 
tion (1+5)  to  give 

CtIo'''  -  .   ^  Cn  -rla  +  ap  -  OiVcos  a  +  .   ^  0.88  sin^(a'  +  acJcos  a' 
■^i  J    A  +  1  "'^  [J        A  +  1  ^      ^^ 

In  view  of  the  very  approximate  nature  of  equation  (1+8)  it  is  desir- 
able to  examine  the  effect  of  this  crossflow  term  on  the  total-lift  coef- 
ficient.  For  the  T\ilin-Burkart,  aspect-ratio-1  section  (A]_  =  0.2) 
used  in  this  investigation,  at  d/c  =  O.O7I  the  ratio  of  the  calcu- 
lated crossflow  lift,   Cl  c  to  the  calculated  total  lift  was  0.157 

at  a  =  1+°  and  O.283  at  a  =  20°.   For  the  five-term  section  of  refer- 
ence 13,  where  Aq^  is  0.075>  the  aspect  ratio  is  3,  and  d/c  is  O.O7I, 
the  ratio  has  been  calculated  as  O.OlU  at  a  =  1+°  and  O.O72  at  a  =  20°. 
Thus  any  inacc\iracies  in  the  crossflow  lift  as  computed  by  equation  (1+8) 
will  appreciably  affect  the  total- lift  coefficient  at  large  angles,  small 
aspect  ratios,  and  large  cambers.   On  the  other  hand  at  higher  aspect 
ratios  and  small  cambers,  errors  in  the  crossflow  component  do  not  greatly 
influence  the  total  lift. 

Equation  (1+9)  may  be  written  in  terms  of  the  slope  m  (given  in 
fig.  2)  as 

o_  _   A   ^/   .  _    _.\     cos  a 


i 


i(cL   +  a^  -  (Xij 


A  +  1  V     ^   ^-^  cos  (a  +  a^  -  Oi) 

0.88  sin^(a'  +ac)cos  a'  (5O) 


A  +  1 


where  ac  is  obtained  from  figure  11  for  the  depth-chord  ratio  of 
interest  and  a^  is  obtained  from  equation  {k6)  .      In  equation  (1+6), 
^L  1  ^^   ^^®  first  term  in  equation  (50).   Equation  (50)  is  solved  by 
iteration  and  the  convergence  is  quite  rapid.   A  sample  calculation  of 
the  lift  coefficient  of  a  cambered  hydrofoil  is  given  in  appendix  B. 
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Drag 

The  drag  coefficient  of  a  supercavitating  hydrofoil  of  finite  aspect 
ratio  operating  at  zero  cavitation  n-omber  and  finite  depth  of  submersion 
is 


^D  ~  '"L  1  ■tan(a  +  e)  +  ^i  ^   tan  a'  +  Cf 


(51) 


where  C^  x  is  the  first  term  in  equation  (50)  and  e   is  the  deviation 
of  the  res'ultant-force  vector  from  the  normal.   For  a  flat  plate  e  =  0, 
a'  =  a,  and  thus  Cq  =  Cl  tan  a  +  Cf .   For  cambered  surfaces  similar  to 
the  circular-arc  or  Tulin-Biirkaxt  section,   €  becomes  very  small  at  large 
angles  of  attack  and  may  be  neglected;  however,  at  small  angles  of  attack, 
the  effect  of  e  on  the  drag  coefficient  cannot  be  neglected.   An  approxi- 
mation to  the  value  of  e  can  be  made  by  determining  its  value  from  the 
two-dimensional  linearized  solution  and  then  modifying  the  resiilt  for  the 
case  of  finite  angles  of  attack  and  aspect  ratio. 

The  linearized  drag  coefficient  as  given  in  reference  1  is 


,^[Z 


Cd=  - 


kxvu    ,- 
cV^ 


(52) 


replacing 
becomes 


X  by 


f(l 


cos  9)   and  die  by 


sin  9  d9  equation  (52) 


Cd  =  - 


(1  -  cos  9) sin  9  uv 


d9 


(53) 


Now  u  may  be  written  in  terms  of  the  vorticity  on  the  equivalent  air- 
foil operating  at  finite  depth  as 


u  =  | 


^0 


1  +  cos  9 
sin  9 


Cn  sin  n9 


n=l 


(5^) 
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and 


V  =  V  -Aq  +  )   A^  cos  ne 

n=l 


(55) 


Therefore 


Cj)  =:  -/   (1  -  cos  9)  sin  9 
^0 


^"(^^fF)^!^"-"-!-*"*! 


-Aq  +  )    ^  COS  n9ld9 
(56) 


For  the  condition  where  A^,  =  a,  C-q     becomes  after  Integrating 


Cd  =  H^o  +  ci  -  ^  U  I 


^  -  ?)Al  ^ 


Co  +  Ci  +  ^ 


2  L    ^^^2 

-yA2  -  -IT 


(57) 


At  infinite  depth   (Cq  =  Aq,      Cn  =  A^)   equation  (57)  reduces  to  the 
value  given  by  Tulin-Burkart  in  reference  1. 

-  \2 


\--i{^*¥)  -s(5J 


(58) 


For  Aq  =  CL,   that  is  Aq'  =  0 


a  +  e  = 


Cd 


f^(^O^C,-^)H.|L(^-^)A, 


Co  +  Cl  +  §t] 


/A2 


A4C2 


k    J 


f  (co  .  c,  -  ^) 


Therefore 


€  = 


(Ci  -  2C2)ai  +  (2C0  +  2Ci  +  Ci^)A2  -  A4C2 


(59) 


(60) 


The  value  of  e  given  by  equation  (60)  is  adequate  only  for  the  case  of 
small  angle  of  attack  and  camber  and  depth-chord  ratios  larger  than 
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about  1.   At  small  spray  thickness  to  chord  ratios  it  may  he  seen  in 
figijre  13  that  the  spray  angle  becomes  quite  large  even  for  small  angles 
of  attack.   When  such  a  flow  is  transformed  by  Z  =  -^zT  the  cavity 
streamline  and  the  free  water  surface  are  rotated  as  shown  in  figure  I3. 
In  figure  13(b)  the  boundary  conditions  are  now  different  from  the  simple 
u  =  0  used  in  the  small-angle  theory.   If  a  system  of  vortices  could  be 
located  to  satisfy  the  boundary  conditions  along  these  new  lines,  a  solu- 
tion for  the  resulting  vorticity  on  the  foil  could  be  obtained.   Such  a 
method  would  involve  taking  a  different  spray  angle  (from  fig.  3)  for 
each  depth  and  angle  of  attack  in  order  to  locate  the  image  vortex  or 
array  of  vortices.   Also,  for  large  angles  of  attack  and  spray  angles 
the  linearizing  assumption  of  v  «  V  is  not  adequate.   In  calculating 
the  lift,  these  difficulties  were  avoided  by  using  Green's  solution  which 
takes  the  effect  of  the  spray  angle  into  account.   It  can  be  seen  in 

Ct 

figure  2  that  as  the  angle  of  attack  increases,  the  ratio  of    ^       dimin- 

ishes,  therefore,  the  image  vortex  must  have  less  influence  on  the 
resulting  hydrofoil  vorticity  for  large  angles  of  attack.   An  approxima- 
tion to  the  correct  hydrofoil-vorticity  distribution  for  finite  angles 
of  attack  operating  near  the  free  water  surface  can  be  obtained  by  using 
the  model  shown  in  figure  8,  and  increasing  the  value  of  d/c  used  so 

Ct 

that  the  resulting    -^   corresponds  to  that  given  by  Green.   Thus,  for 

large  angles  of  attack,  camber,  and  finite  aspect  ratio  an  effective 
depth  of  submersion,   (d/c)g  should  be  used  in  equation  (60) .   The  value 
of   (d/c)g  is  the  value  of  d/c  on  the  a  =  0  line  in  figure  2  corre- 
sponding to  the  value  of  m  =  m^,  where 


me  =  ^  I  =  'f  "  "^  -  "^}   i  (61) 


Cl,oo  2   ni|a  +  a^  -  0^)00  ^ 


The  value  of  the  C  coefficients  are  then  determined  for  (d/c)e  and 
Aq  =  a  -  a^  . 

A  sample  calculation  of  the  drag  coefficient  is  given  in  appendix  B. 
It  has  been  found  after  several  calculations  that  the  value  of  e  and 
thus  the  drag  coefficient  is  not  greatly  affected  by  the  depth  of  sub- 
mersion.  In  fact,  a  rough  approximation  to  e  may  be  obtained  by  assuming 
Cq  =  a  -  cLj_,    and  C^  =  Aj^  in  equation  (60). 
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Center  of  Pressure 

The  linearized  expression  for  the  center  of  pressure  of  a  finite- 
aspect-ratio^  supercavitating  hydrofoil  operating  at  zero  cavitation  num- 
ber and  finite  depth  of  submersion  is 

>'=.p.,l-Ci.i-l6  C2  '^^' 

Cq  +  Ci  -  — 

where  the  C  coefficients  are  determined  at  the  effective  depth  of  sub- 
mersion given  by  eq^oation  (61)  and  for  Aq  =  a  -  o^.   Superimposed  on 
this  flow  is  the  crossflow  component  of  lift  which  is  assumed  to  be 
distributed  uniformly  over  the  chord  and  acting  in  a  direction  normal  to 
the  chord  line.   Thus,  the  distance  from  the  leading  edge  to  the  center 
of  pressure  of  the  crossflow-lift  component  x^^p^  (,      is  given  by 

^c.p.,c  =  0-5c  (63) 

Admittedly,  this  assumption  is  crude  and  accurate  only  for  a  flat  plate. 
For  cambered  surfaces  the  crossflow  will  not  be  uniformly  distributed 
and  for  low  drag  cambered  sections  is  probably  concentrated  on  the  rear- 
ward portion  of  the  hydrofoil. 

By  combining  equations  (62)  and  (63)  the  center  of  pressure  of  the 
combined  flows  is,  therefore, 

/CL,lXc.p.,l  +  0.5Cl,c\  .... 

-c.p.  =  [ cl j'^  (^^ 

A  sample  calculation  of  the  center  of  pressure  is  given  in  appendix  B. 
As  in  the  case  of  e,  a  few  calculations  reveal  that  a  fair  approximation 
for  Xq_p_  ]_  is  obtained  by  using  Cq  =  a  -  aj_  and  C-^   =  A^^  in  equa- 
tion (62)  . 


EXPERIMENTAL  INVESTIGATION 


MODEI£ 


Two  models  of  7-071-inch  chord  were  used  in  the  experimental  inves- 
tigation. As  shown  in  figure  li+(a),  the  first  had  a  lower  surface 
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conforming  to  the  T-ulin-Burkart  low-drag  configuration  given  by  equa- 
tion (12) .   The  two-dimensional  design-lift  coefficient  was  selected  as 
0.392  corresponding  to  a  value  of  0.2  for  the  coefficient  A]_.   Since  the 
hydrofoil  is  designed  to  operate  in  a  cavity,  the  shape  of  the  upper  sior- 
face  is  arbitrary  as  long  as  it  does  not  interfere  with  the  formation  of 
the  cavity  from  the  leading  edge.   Since  the  greatest  advantage  is  to  be 
obtained  at  small  angles  of  attack  ajid  thus  thin  cavities,  the  thickness 
of  the  hydrofoil  must  be  small.   For  the  present  cambered  model  the  upper 
surface  profile  from  leading  edge  to  midchord  was  arbitrarily  chosen  to 
conform  with  the  free  streamline  leaving  the  leading  edge  of  a  two- 
dimensional  flat  plate  at  5°  incidence.   (See  ref.  1.)   The  thickness  of 
the  portion  rearward  of  the  midchord  was  made  in  the  image  of  the  forward 
portion  and  resulted  in  a  symmetrical  thickness  distribution  with  a  maxi- 
miom  thickness-chord  ratio  of  J- 3  percent.   Since  the  center  of  pressure 
of  the  Tulin-Burkart  hydrofoil  is  located  near  the  midchord,  this  symmet- 
rical section  minimizes  the  torsional  moment  on  the  foil  and  results  in 
less  twist  than  wo\ild  be  experienced  by  a  nonsymmetrical  section. 

The  second  model  was  of  triangular  cross  section  with  a  flat  bottom 
as  shown  in  figure  14(b) .   The  maximum  thickness  was  5  percent  of  the 
chord. 

Such  thin  sections  lead  to  structural  limitations  in  the  aspect  ratio 
when  supported  by  a  single  strut  at  midspan.   Since  an  aspect  ratio  of  1 
is  about  the  most  desirable  from  the  structixral  standpoint  and  it  also 
represents  the  accepted  dividing  line  between  hydrofoils  and  hydroskis, 
both  models  were  made  with  a  square  plaji  form. 

The  strut,  which  can  also  be  seen  in  figure  Ik,   had  an  NACA  66]_-012 
airfoil  section.   The  strut  was  mounted  perpendic\ilar  to  the  flat  plate 
and  perpendicular  to  the  X-axis  of  the  cambered  surface.   The  intersection 
of  the  strut  and  upper  surface  was  without  fillets.   Both  the  hydrofoil 
and  the  strut  were  made  of  stainless  steel  and  were  polished  to  a  smooth 
finish. 


APPARATUS  AND  PROCEDURE 


Tests  were  made  by  using  the  carriage  and  existing  strain-gage  bal- 
ances in  the  Langley  tank  no.  2  which  independently  m.easure  the  lift, 
drag,  and  pitching  moment.   Fig\ure  I5  shows  a  view  of  the  test  setup  with 
the  cambered  hydrofoil  and  the  balance  attached  to  the  carriage.   The 
moment  was  measured  about  an  arbitrary  point  above  the  model  and  the  data 
thus  obtained  were  used  to  calculate  the  moments  about  the  leading  edge. 
The  positive  directions  of  forces,  angles,  and  moments  used  in  presenting 
the  force  data  are  shown  in  figure  16. 
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The  force  and  moment  measurements  were  made  at  constant  speeds  for 
fixed  angles  of  attack  and  depths  of  submersion.   The  depth  of  submersion 
is  defined  as  the  distance  from  the  undisturbed  water  surface  to  the 
leading  edge  of  the  model. 

Ventilated  or  zero-cavitation-number  tests  were  made  on  both  models 
at  a  depth  of  submersion  of  O.5  inch  throughout  a  range  of  angle  of  attack 
from  6°  to  20°  for  the  flat  plate,  and  8°  to  20°  for  the  cambered  hydro- 
foil.  Three  methods  of  obtaining  ventilated  flow  at  this  0.5-inch  depth 
were  used:   (l)  Normal  ventilation  through  the  trailing  vortices  as 
described  in  reference  h,    (2)  injection  of  air  through  the  port  on  the 
strut  leading  edge  (see  fig.  li|)  (this  air  was  supplied  at  a  rate  of 
0.012  pounds  per  second  and  cut  off  after  ventilation  was  established), 
and  (3)  a  l/32-inch-diameter  wire  was  soldered  on  the  leading  edge  of 
the  cambered  model  to  cause  local  separation  and  thus  ventilation.   Both 
models  were  also  investigated  at  a  depth  of  submersion  of  0  inch  at 
a  =  1|°  for  the  flat  plate,  and  a  =  6°,  8°,  and  10°  for  the  cambered 
surface.   At  angles  of  attack  of  l6°  and  20°  for  the  flat  surface,  forces 
were  also  measured  over  a  range  of  depth  of  submersion  for  which  venti- 
lation could  be  obtained  (d  =  0  to  2  inches) . 

The  thickness  and  direction  of  the  jet  or  spray  leaving  the  leading 
edge  of  the  flat  plate  were  also  measured  at  l6°  and  20°  for  a  range  of 
d  from  0  to  2  inches.   A  schematic  drawing  of  the  instrument  used  for 
measuring  the  spray  thickness  and  direction  is  shown  in  figure  I7.   The 
stagnation  tube  was  lowered  through  the  spray  during  a  test  and  the  pres- 
sure and  location  of  the  tube  center  line  recorded  on  an  oscillograph. 
Almost  instantaneous  response  of  the  st agnation -tube — pressure-cell 
combination  was  obtained  by  completely  filling  the  tube  and  connecting 
line  with  water.   The  point  of  entering  and  leaving  the  spray  was  obtained 
by  comparing  the  location  of  the  tube  with  the  rise  and  fall  of  pressure 
as  the  tube  passed  through  the  spray.   The  vertical  location  of  the  tube 
was  obtained  from  the  output  of  the  slide -wire  circuit  also  shown  in  fig- 
ure 8.   During  a  test  the  tube  was  passed  through  the  spray  several  times 
and  an  average  of  the  results  was  taken. 

Tests  of  the  nonzero-cavitation-number  case  were  made  at  a  depth  of 
submersion  of  6  inches  where  vortex  ventilation  did  not  occur.   Measirre- 
ments  of  lift,  drag,  and  moment  were  obtained  for  a  range  of  velocities 
from  20  to  80  feet  per  second  at  angles  from  l6°  to  20°  for  which  long 
trailing  air  or  vapor  cavities  could  be  obtained.   Data  for  air-filled 
cavities  were  obtained  by  introducing  air  from  an  external  metered  supply 
to  the  upper  surface  through  the  ports  on  the  strut  leading  edge  shown 
in  figure  lU.   During  these  tests  the  pressure  in  the  cavity  formed  was 
measured  by  a  pressure  cell  connected  to  a  l/l6-inch-diameter  orifice 
located  near  the  bottom  of  the  strut  so  as  to  be  within  the  cavity.   (See 
fig.  lii.)   This  measured  pressure  was  used  in  computing  the  cavitation 
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number  for  the  cavity  formed.   The  airflow  rate  was  measured  by  an 
orifice-type  flow  meter.   Airflow  rates  up  to  0.012  poimd  per  second 
were  obtained  with  the  test  arrangement. 


ACCURACY 


The  change  in  angle  of  attack  due  to  structural  deflection  caused 
by  the  forces  on  the  model  was  obtained  during  the  calibration  of  the 
balances  and  test  data  were  adjusted  accordingly.   The  maximum  correction 
necessary  was  only  0.1  .   The  estimated  accuracy  of  the  measurements  is 
as  follows : 

Angle  of  attack,  deg ±0.1 

Depth  of  submersion,  in ±0.1 

Speed,  fps ±0.2 

Lift,  lb ±15 

Drag,  lb ±7 

Moment,  ft-lb   ±6 

Cavity  pressure,  Ib/sq  ft ±10 

Spray  thickness,  in ±0.05 

Spray  angle,  deg ±1.5 

The  forces  and  moments  were  converted  to  the  usual  aerodynamic  coef- 
ficient form  by  using  a  measured  value  of  the  density  of  1.93  slugs  per 
cubic  foot.   The  kinematic  viscosity  measured  during  the  tests  was 

1.70  X  10   pounds-second  per  square  foot.   Thus,  for  the  range  of  veloc- 
ities investigated,  the  Reynolds  number  based  on  chord  ranged  from 

0.7  X  10^  to  2.8  X  10^. 


RESULTS  AND  DISCUSSION 


ZERO  CAVITATION  NUMBER 


General 

Ventilation  inception.-  The  process  of  ventilation  by  air  entrainment 
through  the  tip  vortices  of  a  hydrofoil  near  the  free  water  siirface  is 
described  in  reference  k   and  is  shown  for  the  present  models  operating  at 
a  depth  of  O.5  inch,  in  figures  I8  and  I9.   The  speed  at  which  complete 
ventilation  occurred,  when  tested  at  this  depth  is  shown  in  figiire  20. 
Complete  ventilation  by  entrainment  through  the  vortices  was  not  possible 
at  angles  less  than  10°  for  the  flat  surface  and  12°  for  the  cambered 
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surface.   Figure  l8(a)  shows  the  incomplete  ventilation  of  the  flat  plate 
at  8°  and  figure  19(a)  the  similar  action  of  the  ca^-faered  surface  at  10  . 
Figure  iBCb)  shows  the  complete  ventilation  of  the  flat  surface  at  10°  and 
figures  19Cb)  and  (c),  the  complete  ventilation  of  the  cambered  surface 
at  12°  and  16°,  respectively.   Leading-edge  cavitation  of  hoth  models  may 
"be  noted  in  the  photographs. 

One  of  the  requirements  for  establishing  ventilated  flow  from  the 
leading  edge  is  that  the  hydrofoil  upper  surface  lie  beneath  the  upper 
cavity  streamline  for  the  condition  of  zero  cavitation  number.   Thus, 
for  a  given  hydrofoil  thickness  there  is  a  minimum  angle  below  which 
ventilation  from  the  leading  edge  caimot  exist.   Figure  21  shows  the 
theoretical  infinite  depth  location  of  the  upper  cavity  boundary  for  two- 
dimensional  hydrofoils  having  the  same  cross  section  as  the  models  tested. 
These  streamlines,  computed  from  equations  given  in  references  1  and  I5, 
reveal  that  ventilation  of  the  two-dimensional  flat  plate  should  be 
possible  at  angles  greater  than  3.2  . 

As  pointed  out  in  references  1  and  I5,  the  slope  of  the  cavity 
streamline  is  determined  by  the  magnitude  and  distribution  of  the  circu- 
lation on  the  equivalent  airfoil  section.   The  reduction  of  the  two- 
dimensional  (infinite  depth)  circulation  by  an  amount  equivalent  to 
Jones'  edge  correction,   1/E  =  A/(A  +  l) ,  will  cause  the  cavity  streamline 
to  be  lowered  over  the  entire  span  by  an  amount  proportional  to  l/t.. 
When  the  induced  angle  oi  Is  considered,  the  streamline  is  brought  even 
closer  to  the  hydrofoil  upper  surface,  particularly  near  the  tips. 

By  considering  these  two  effects  and  by  assuming  that  the  influence 
of  oi  is  uniform  over  the  span,  the  effect  of  finite  aspect  ratio  on 

^^\       may  be  approximated  by  the  following  equation: 


the  cavity  slope 


■dy\ 


■dy' 
dxVc 


(,a  -  ai) 

A  +  1 

(a')A=oo 


A=oo 


(65) 


Thus,  in  order  to  establish  a  cavity  at  angle  of  attack  a  and  aspect 
ratio  A  which  is  equivalent  to  the  cavity  formed  at  angle  of  attack  a 
and  A  =  00  the  following  equation  must  be  satisfied: 


a 


~  A  + 


-(a-)  A., 


A=oo 


+  aj_ 


(66) 
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By  using  equation  (66)  and  the  value  of   (a-);^_^  =  3.2°,  the  geometric 
angle  of  attack  of  the  aspect-ratio-1  plate  would  have  to  be  about  8° 
in  order  to  ventilate  from  the  leading  edge  at  infinite  depth  of  sub- 
mersion.  The  effect  of  the  proximity  of  the  free  water  s\irface  is  to 
rotate  the  cavity  streamline  away  from  the  hydrofoil  as  seen  in  fig- 
ure 3.   Thus,  it  should  be  possible  to  establish  a  vented  cavity  on  the 
flat  section  at  angles  at  least  as  small  as  8°  when  operating  near  the 
free  water  surface.   Therefore,  the  absence  of  ventilation  at  angles 
less  than  10°  cannot  be  attributed  to  interference  from  the  model  upper 
surface.   In  fact  the  formation  of  leading-edge  cavitation  at  8°  indi- 
cates that  if  the  velocity  could  have  been  increased  well  above  80  feet 
per  second  (preferably  to  zero  cavitation  number)  a  vapor  cavity  com- 
pletely inclosing  the  upper  surface  would  result.   The  reason  that  venti- 
lation did  not  occur  at  angles  less  than  10°  could  be  that  the  region  of 
boundary- layer  separation  in  the  midportion  of  the  chord  did  not  extend 
to  the  aerated  vortices.   This  conjecture  was  proved  to  be  the  case  by 
injecting  air  from  the  external  supply  down  the  strut  to  the  upper  sur- 
face of  the  model.   Complete  ventilation  of  the  upper  surface  was  then 
possible  at  angles  as  low  as  6°.   The  air  was  supplied  (at  the  maximum 
rate  of  0.012  pound  per  second)  during  the  acceleration  of  the  model 
and  shut  off  after  ventilation  was  established.   The  process  is  shown 
in  figure  l8(c)  for  the  flat  plate  at  8°  angle  of  attack. 

Figure  21  shows  that  ventilation  of  the  two-dimensional  cambered 
section  should  be  possible  at  angles  above  about  4°.   In  reference  I3 
at  large  angles  of  attack,  the  cavity  shape  is  almost  independent  of  the 
camber,  depending  only  on  the  angle  of  attack.   Therefore,  the  finite- 
aspect-ratio  corrections  given  by  equations  (65)  and  (66)    are  assumed 
applicable  to  cambered  sections.   When  these  corrections  for  finite 
aspect  ratio  are  taken  into  account  the  minimum  angle  for  ventilation  of 
the  cambered  section  becomes  about  12°.   Again  the  proximity  of  the  free 
water  surface  will  reduce  this  infinite-depth  estimate.   Thus,  failure 
to  establish  complete  ventilation  of  this  hydrofoil  at  angles  less  than 
12°  (at  d  =  0.5  inch)  may  have  been  due  to  interference  of  the  upper 
boundary  of  the  model.   In  fact  it  was  found  that  injecting  air  to  i^he 
upper  surface  at  angles  less  than  12°  failed  to  establish  ventilation 
or  to  influence  the  ventilation  speed  at  12°.   Since  the  lift-drag  ratio 
of  cambered  supercavitating  hydrofoils  increases  with  decrease  in  angle 
of  attack,  these  observations  suggest  the  use  of  higher  geometric  aspect 
ratios  or  the  effective  increase  of  aspect  ratio  by  the  use  of  end  plates. 
End  plates  appear  to  be  particularly  applicable  if  low-aspect-ratio 
cambered  hydrofoils  are  to  be  designed  to  ventilate  from  the  leading  edge 
at  small  angles  of  attack  and  large  depths  of  submersion. 

Since  complete  ventilation  seemed  to  require  extensive  boundary- 
layer  separation  on  the  upper  surface,  a  l/32-inch  diameter  wire  was 
soldered  on  the  leading  edge  of  the  cambered  model  to  cause  a  local 
boundary- layer  separation.  With  this  wire  in  place  it  was  possible  at 
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the  0.5-inch  depth  of  submersion  to  obtain  complete  vortex  ventilation 
of  the  upper  surface  at  10°  incidence  but  not  at  8  incidence  even  at  the 
maximum  available  velocity  of  80  feet  per  second.   However,  by  introducing 
air  through  the  ports  on  the  strut  it  was  possible  to  obtain  complete 
ventilation  at  angles  as  low  as  8°.   Evidently  air  brought  up  by  the  tip 
vortices  at  8°  incidence  could  not  reach  the  small  separated  region  behind 
the  wire;  however,  the  forced  air  was  capable  of  reaching  this  zone  and 
establishing  a  cavity.   The  cavity  formed  extended  downstream  far  enough 
to  intercept  the  ventilated  tip  vortices  and  complete  ventilation  then 
occurred.   Of  course,   the  primary  influence  of  the  wire  was  to  increase 
the  cavity  ordinates  so  that  the  hydrofoil  upper  surface  did  not  inter- 
fere with  establishing  ventilation  from  the  wire.   Although  the  basic 
force  data  obtained  with  the  wire  were  recorded  and  are  presented,  they 
are  not  used  in  the  later  correlation  with  theory  because  of  the  unknown 
influence  of  the  wire  on  the  forces . 

Basic  force  and  moment  results . -  The  basic  data  from  the  tests  of 
the  ventilated  hydrofoils  at  a  depth  of  submersion  of  O.5  inch  are  pre- 
sented in  figure  22  for  the  flat  plate  and  figure  23  for  the  cambered 
s\irface  as  curves  of  lift,  drag,  and  pitching  moment  about  the  leading 
edge  against  speed,  for  various  angles  of  attack.   Ventilated-flow  data 
obtained  at  zero  depth  of  submersion  at  k'-'   incidence  for  the  flat  plate 
and  6°,  8°,  and  10°  for  the  cambered  surface  are  also  included. 

The  basic  ventilated  flow  data  obtained  for  depths  of  0,  1.0,  1.^, 
and  2.0  inches  for  incidences  of  l6°  and  20°  &re   presented  in  figure  2k, 


Comparison  of  Experimental  Results  With  Theory 

Spray  thickness. -  In  the  section  on  theory  the  need  for  determining 
the  relationship  between  the  leading-edge  depth  of  submersion  and  the 
spray  thickness  for  a  flat  plate  was  pointed  out.   These  variables  were 
measured  for  the  flat  plate  at  l6°  and  20°  over  the  range  of  depth  of 
submersion  from  0  to  2  inches.   The  data  are  shown  in  figure  25.   The 
spray-thickness  measurements  presented  were  obtained  approximately  20  per- 
cent of  the  chord  rearward  of  the  leading  edge.   On  one  test,  measurement 
of  the  spray  thickness  was  also  obtained  at  about  the  midchord  and  the 
results  were  in  agreement  with  those  obtained  at  the  20-percent  location. 
Therefore,  for  the  range  of  depths  investigated,  the  spray  thickness 
measured  can  be  considered  as  the  theoretical  value  infinitely  rearward 
of  the  foil.   It  may  be  noted  in  figure  25  that  the  spray  thickness  is 
greater  than  the  leading-edge  depth  of  submersion,  the  magnitude  of  the 
ratio  increases  with  either  increase  in  angle  of  attack  or  decrease  in 
depth  of  submersion.   Also  shown  in  figure  25  are  the  theoretical  rela- 
tionships between  6/c  and  d/c  for  the  two-dimensional  0°  and  90°  cases 
previoiisly  discussed.   It  may  be  noted  that  the  trends  of  the  experimental 
and  theoretical  curves  are  parallel. 
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In  foregoing  sections  it  has  been  pointed  out  that  the  relationship 
between  leading- edge  depth  of  submersion  and  spray  thickness  is  influ- 
enced by  the  force  of  gravity  and  thus  cannot  be  obtained  from  Green's 
analysis.   Also  the  influence  of  aspect  ratio  has  been  shown  to  be  impor- 
tant.  The  relationship  for  a  =  0°  and  90°  has  been  determined  but  the 
manner  in  which  intermediate  angles  of  attack  affect  the  result  was  not 
understood.   The  data  shown  in  figure  25  combined  with  the  end-point 
results  previously  established  for  a  =  0°  and  90°  permit  lines  of  con- 
stant d/c,  at  least  for  the  aspect-ratio-1  condition,  to  be  drawn  on 
figure  2.   The  experimental  data  shown  in  figure  25  were  plotted  on  fig- 
ure 2  at  their  equivalent  angle  of  attack  a  -  a^ .   This  corresponds  to 
equivalent  angles  of  attack  of  approximately  12.3°  fo^  the  l6°  case  and 
15.9°  for  the  20°  case  for  the  range  of  d/c  presented.   For  values  of 
d/c  greater  than  O.285,  only  the  theoretical  end  points  were  available 
and  paralleling  lines  were  faired  in.   Although  these  lines  are  accurate 
only  for  the  aspect-ratio-1  condition,  they  are  considered  to  be  good 
approximations  even  for  aspect  ratios  as  high  as  6.   Such  an  approxima- 
tion is  reasonable  because  the  asymptotic  value  of  the  stagnation  line 
infinitely  forward  of  a  flat  plate  is  only  of  the  order  of  about  0.1  chord 
below  the  stagnation  point  for  aspect  ratios  less  than  6  and  lift  coef- 
ficients less  than  0.5-   It  is  obvious  from  figure  2  that  for  depths 
greater  than  1  chord,  the  end  points  may  be  connected  by  any  reasonable 
line  (for  example,  a  straight  line)  with  very  little  loss  in  accuracy. 

Spray  angle .  -  Figiore  26  shows  the  effect  of  depth  of  submersion  on 
the  spray  angle  (f)     for  the  flat  plate  at  l6°  and  20°.   The  theoretical 
spray-angle  variation  also  shown  in  figure  26  was  obtained  from  figure  5 
by  considering  the  two-dimensional  angle  of  attack  as  a  -  ai  and  the 
two-dimensional  6/c  as  the  actual  measured  value.   The  tangent  of  the 
spray  angles  obtained  from  figure  3  for  the  angle  a  -  cuj^  were  then 
reduced  by  the  factor  1/E.   The  justification  for  this  1/E  modification 
is  the  same  as  that  presented  in  the  section  on  ventilation  inception. 
That  is,  that  the  slope  of  the  upper  caA/lty  streamline  is  proportional 
to  the  equivalent  airfoil  circulation  and  thus  if  the  two-dimensional 
circulation  is  reduced  by  l/E,  the  two-dimensional  spray-angle  slope 
will  also  be  reduced  by  this  amount.   When  these  three-dimensional  cor- 
rections are  made,  the  calculated  angles  are  still  high  at  small  values 
of  d/c  and  appear  to  coincide  with  the  experimental  data  for  values 
of  d/c  greater  than  about  O.25.   Since  the  theorem  that  the  slope  of 
the  cavity  streamline  is  directly  proportional  to  the  circulation  was 
obtained  in  reference  1  for  very  thin  cavities  at  infinite  depth;  it 
cannot  be  expected  to  be  applicable  when  the  spray  angles  are  very  large. 
This  deviation  from  assumed  small  angles  may  be  the  only  explanation 
req;iired  for  the  discrepancy  between  theory  and  experiment  at  small  values 
of  d/c.   However,  it  might  also  be  attributed  to  the  influence  of  gravity 
and  finite  aspect  ratio  on  Green's  solution.   In  the  section  entitled 
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"Relation  Between  Depth  of  Submersion  and  Spray  Thickness,"  it  was  sug- 
gested that  the  spray  angles  given  in  figure  3  would  be  reduced  because 
of  these  effects  and  at  least  this  assumption  has  not  been  disproved  by 
the  experimental  data. 

Lift  coefficient.-  All  ventilated  force  and  moment  data  in  coefficient 
form  were  found  to  be  independent  of  speed  in  the  range  tested.   (This 
independence  was  true  in  the  present  investigation  because  of  the  shallow 
depth  of  submersion  and  therefore   0=0.   At  very  large  depths,  a     will 
be  greater  than  zero  because  Pq  is  greater  than  p^  even  if  the  cavity 
is  fully  vented  to  the  atmosphere.   Therefore  at  large  depths,  changes  in 
velocity  will  affect  the  lift  coefficient  because  these  changes  affect  a. ) 
The  data  shown  in  figures  22(a)  and  23(a)  are  plotted  in  figure  27  as  lift 
coefficient  against  angle  of  attack  for  each  of  the  models  tested.   A 
comparison  of  the  lift  coefficients  of  the  two  models  shows  an  effective 
increase  in  angle  of  attack  of  the  cambered  model  as  predicted  by  the 
Tulin-Burkart  theory.   Also  shown  in  figure  27  are  the  theoretical  lift- 
coefficient  ciurves  obtained  from  equation  (50)-   The  theory  is  about 
5  percent  lower  than  the  measured  values.   Thus,  the  use  of  equation  (50) 
in  engineering  calculations  of  the  lift  coefficient  appears  to  be  warranted. 

The  variation  of  lift  coefficient  with  depth  of  submersion  of  the 
flat-plate  model  at  angles  of  incidence  of  l6°  and  20°  is  shown  in  fig- 
ure 28.  Note  the  slight  increase  in  lift  coefficient  as  the  hydrofoil 
nears  the  surface.  Also  shown  in  figure  28  is  the  theoretical  variation 
of  the  lift  with  depth  of  submersion  obtained  from  equation  (50)*  The 
theory  is  in  excellent  agreement  with  the  data  and  accurately  predicts 
the  increase  in  lift  with  decreasing  depth  of  submersion. 

Further  verification  of  equation  (50)  is  shown  in  figure  29  where  a 
comparison  is  made  with  the  experimental  data  of  Fuller  (ref .  I7) .   These 
data  were  obtained  on  sections  with  lower  surface  profiles  similar  to 
those  of  the  present  investigation  but  at  an  aspect  ratio  of  2.   The 
calculated  values  are  in  good  agreement  with  the  measured  values  for 
both  the  flat  and  cambered  models . 

Drag  coefficient.-  In  figure  50  the  data  of  figures  22(b)  and  23(b) 
in  coefficient  form  are  compaxed  with  theoretical  values  obtained  from 
equation  (51)-   The  friction  drag  coefficient  of  one  side  of  either  of 
the  models  was  calculated  to  be  about  O.OO3.  With  the  strut  drag  included, 
the  total-skin-friction  drag  coefficient  Cf     was  taken  as  O.OO^i.   By 
using  this  value  of  Cf  in  equation  {^l) ,   the  agreement  between  theory 
and  experiment  is  good  for  both  models. 

In  figure  3I  the  experimental  lift-drag  ratios  obtained  from  the 
data  of  figures  22  and  23  are  compared  with  theory.   Again,  both  experi- 
ment and  theory  include  the  skin-friction  drag  coefficient   (Cf  =  O.OOU). 
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The  agreement  between  theory  and  experiment  is  good  for  both  models. 
The  superiority  of  the  cambered  hydrofoil  is  clearly  revealed  in  this 
figure.   At  a  lift  coefficient  of  0.25  the  L/d  of  the  cambered  hydro- 
foil is  more  than  twice  that  of  the  flat  plate. 

Also  included  in  figure  3I  sJ^e  data  taken  from  reference  if  on  a 
ventilated  modified  flat  plate  of  aspect  ratio  1.   This  modified  plate 
had  an  elliptical  nose  and  a  tapered  trailing  edge.   The  importance  of 
providing  a  sharp  leading  edge  on  hydrofoils  designed  for  use  in  cavity 
flow  is  shown  by  comparing  the  l/D  of  this  modified  flat  plate  with 
the  L/d  of  the  sharp-nosed  flat  plate  of  the  present  investigation. 
The  rounded  leading  edge  of  the  modified  plate  is  subjected  to  a  net 
positive  pressure  which  is  not  balanced  by  similar  pressures  on  the  rear- 
ward portion  of  the  plate.   In  addition,  the  lower  surface  of  this  plate 
presents  an  effective  negative  camber  to  the  flow  axid  thus  does  not  develop 
as  much  lift  as  a  truly  flat  surface.   At  small  angles  the  drag  of  the 
rounded  nose  greatly  influences  the  maximum  L/D  of  the  section.   Several 
other  investigations  have  noted  the  importance  of  a  sharp  leading  edge  on 
hydrofoils  designed  for  operation  in  the  supercavitating  regime.   (See 
refs.  18  and  I9.) 

Center  of  pressi^re.  -  The  center  of  pressiire  of  the  flat  and  cambered 
models  calculated  from  the  data  of  figiires  22  and  25  are  compared  with 
theory  in  figure  32.   The  theory  from  equation  {6k)    is  in  good  agreement 
with  the  experimental  data  for  both  models.   Since  the  accuracy  of  the 
forces  and  moments  on  the  flat  plate  is  poor  at  small  angles  of  attack 
(small  total  loads)  the  accuracy  of  the  center  of  pressure  from  the  data 
obtained  on  the  flat  plate  at  angles  of  incidence  of  J|°  and  6°  is  doubtful. 


NONZERO  CAVITATION  NUMBER 


The  nonzero-cavitation-number  characteristics  of  the  two  models 
obtained  at  a  6-inch  depth  of  submersion  are  shown  in  figure  33  for  angles 
of  attack  of  l6°,  l8°,  and  20°.   The  solid  data  points  are  the  lift  coef- 
ficients obtained  for  vapor  cavitation.   The  cavitation  number  corresponding 
to  the  condition  tested  was  computed  by  using  the  water  vapor  pressure  at 
the  test  temperature  for  the  pressure  within  the  cavity  /that  is 

0  =  —2 V  ).   For  cavitation  numbers  less  than  about  O.7  the  vapor  pres- 

q   / 
sure  was  the  same  as  the  measured  cavity  pressure.   However,  at  cavitation 

numbers  greater  than  0.7^  the  measured  pressure  was  usually  higher  than 
the  vapor  pressure.   Since  it  could  not  be  determined  with  certainty 
whether  the  cavity  pressure  orifice  was  within  the  cavity,  the  vapor  pres- 
sure was  used  to  compute  all  cavitation  numbers  denoted  by  the  solid  data 
points.   Also  denoted  in  figure  35  are  the  approximate  values  of  a     at 
inception  a^     and  the  point  at  which  the  cavity  length  exceeded  the  chord. 
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These  values  are  only  estimates,  since  no  effort  was  made  to  find  the 
exact  velocity  at  which  these  incidents  occurred.   Although  the  estimated 
values  of  a±     are  the  same  for  both  models,  there  may  act\xally  be  some 
difference  in  the  true  points. 

The  blank  data  points  in  fig\:ire  35  represent  data  obtained  by  intro- 
ducing air  to  the  upper  surface  of  the  model  and  establishing  a  cavity. 
The  cavitation  number  for  this  condition  was  computed  by  using  the 
measured  value  of  the  pressure  within  the  cavity.   In  this  case  the  cavity 
pressure  orifice  was  always  well  within  the  cavity  formed. 

In  figure  35  the  curves  are  drawn  through  the  computed  value  of  C^ 
for  CT  =  0  as  obtained  from  equation  (50) . 

The  agreement  between  the  vapor-  and  air-cavity  data  confirms  the 
use  of  the  cavitation  number  as  the  significant  parameter  for  correlating 
the  characteristics  of  cavity  flow.   The  similarity  of  the  air  and  vapor 
cavities  at  nearly  equal  cavitation  number  is  shown  in  figure  3^-   This 
dependence  on  the  cavitation  number  is  to  be  expected  because  the  forces 
on  the  body  are  influenced  only  by  the  streamline  curvatures  and  thus  the 
pressure  within  the  cavity;  the  type  of  gas  present  should  have  only  a 
secondary  Influence. 

It  was  not  possible  to  establish  satisfactory  air  cavities  at  angles 
less  than  about  l6°.   At  low  speeds  where  the  air  could  reach  the  sepa- 
rated region  near  the  leading  edge  the  cavity  upper  surface  was  greatly 
disturbed  by  the  force  of  the  air  jets.   At  higher  velocities  either  the 
air  could  not  reach  far  enough  upstream  to  form  a  cavity  or  when  a  cavity 
was  formed  it  did  not  cover  the  whole  chord.   If  a  greater  quantity  of 
air  is  supplied  it  is  believed  that  satisfactory  results  can  be  obtained 
at  angles  less  than  l6°. 

The  use  of  a  dynamic  model  of  a  high-speed  aircraft  will  require 
simultaneous  reproduction  of  both  the  full-scale  cavitation  and  Froude 
nijmbers.   Such  an  investigation  is  possible  if  an  air-filled  cavity  with 
the  proper  cavitation  number  can  be  established  on  the  model. 


CONCLUSIONS 


Conclusions  based  on  the  resiilts  obtained  from  the  theoretical  and 
experimental  investigations  of  supercavitating  flat  and  cambered  hydro- 
foils may  be  summarized  as  follows: 

1.  The  theoretical  expressions  derived  for  the  lift,  drag,  and 
center  of  pressure  of  supercavitating  hydrofoils  of  arbitrary  section 
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operating  at  zero  cavitation  niamber,  and  finite  aspect  ratio  and  depth 

of  submersion  are  in  good  enough  agreement  with  the  available  experimental 

data  to  warrant  their  use  in  engineering  calculations. 

2.  The  experimental  ventilated  force  and  moment  coefficients  of  both 
models  investigated  at  shallow  depths  of  submersion  were  independent  of 
speed. 

3.  Similar  to  the  influence  of  camber  on  a  fully  wetted  airfoil,  the 
influence  of  camber  in  ventilated  flow  was  to  effectively  increase  the 
angle  of  attack  without  appreciably  changing  the  lift-curve  slope. 

k.    The  lift-drag  ratio  of  the  cambered  model  near  the  design  lift 
coefficient  was  more  than  twice  that  of  the  flat  plate. 

5.  Comparison  of  the  sharp-nosed  flat  plate  with  a  rounded-nose  flat 
plate  showed  the  sharp-nosed  section  to  be  considerably  superior.   Thus, 
hydrofoils  designed  for  operation  in  the  supercavitating  regime  of  flow 
should  have  sharp  leading  edges  for  best  efficiency. 

6.  The  cavitation  number  defines  the  flow  similarity  and  lift  of  a 
supercavitating  hydrofoil  regardless  of  the  type  of  gas  (air  or  water 
vapor)  within  the  cavity.   Thus  vapor  cavities,  normally  obtained  at  high 
speeds  can  be  simulated  at  lower  speeds  by  establishing  a  cavity  with 
forced  air  whose  cavitation  number  is  the  same  as  that  of  the  high-speed 
flow. 


Langley  Aeronautical  Laboratory, 

National  Advisory  Committee  for  Aeronautics, 
Langley  Field,  Va.,  August  27,  1957- 
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APPENDIX  A 

THE  TULIN-BURKART  LINEAR  THEORY  FOR  ZERO-CAVITATION-NUMBER 
.  FLOW  APPLIED  TO  A  CIRCULAR-ARC  HYDROFOIL 


The  problem  is  to  find  the  effective  increase  in  angle  of  attack  due 
to  the  camber  of  a  two-dimensional  circular-arc  hydrofoil  operating  at 
infinite  depth  and  zero  cavitation  number  by  the  method  of  reference  1, 
that  is,  to  find  the  section  coefficients  Aq,      A2_,    and  A2  in  the  fol- 
lowing equation  if  a,  measured  from  the  chord  line,  is  zero: 


CL  =  f(Ao'  ^Ai-f)  =  f(Ao'  -.aj 


By  using  the  system  of  axes  shown  in  the  following  sketch 
Y 


c 

7 

\  '' 

/ 

c/2 

/ 

and  the  notation  of  reference  1  (barred  symbols  for  the  equivalent  air- 
foil and  unbaxred  for  the  hydrofoil) ,  the  equation  of  the  circiilar-arc 
hydrofoil  is 


(x  -  c/2)^  +  (y  +  R  cos  7/2)'^  =  R 


2  _  T,2 


(Al) 


dy 
dx 


X  -  c/2 
/r2  -  (x  -  c/2) 2 


(A2) 
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From  reference  1  the  slope  of  the  eq\ilvalent  airfoil  is  obtained  from 
the  equation 


^(x)  =  ^{^^)  (A5) 

dx     dx 


_2 

By  noting  that  c  =  c  ,  the  slope  of  the  eqxiivalent  airfoil  section  is, 

therefore, 


dy  ^      x^  -  5^/2 

d5     n;   7T^   72/012 


Jr2  -  (x2  -  52/2)= 


If  the  conventional  substitution  of 


aj?  _     |(i  -  cos  9)  -  1 


dx 


\ 


^-   [|(i-cos  ef  -i 


but 


Therefore, 


r2   r2   /_^\2    /^  v2 
Z'c£"W2y    (sinZ 


2 

dy  _   1       cos  9-2  cos  9-1 

dx     2  /        ^    ~~  r^ 

1     /eos29  -  2  cos  8-1] 

sin2  Z   V        2        ; 


(A4) 


X  =  |(1  -  cos  9)  (A5) 


is  made,  equation  (h)    after  some  manipulation  may  be  written  as 


(A6) 


(AT) 


(A8) 
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Since  the  linear  theory  is  applicable  only  to  small  7  and, 
therefore, 

1    »  /co3^9  -  2  cos  9  -  l\ 


•  2  r 
2 


equation   (A8)    may  "be  approximated  as 


^  ^  _  Zfcos^e   -  2  cos   9   -  1)  (A9) 

dx  l^. 


The  necessary  coefficients  are  then  readily  obtainable  as 


Ai=ir^ 


^  cos  9  d9  =  Z  (All) 

dx  2 


rt  ,- 


A^  =  ^  /   ^  cos  29  d9  =  -  ^  (A12) 


An  =  0  (n  >  2)     (AI5) 


The  required  effective  increase  in  angle  of  attack  due  to  small  circular - 
arc  camber  is,  therefore. 


It  shoiild  be  noted  that  for  the  reference  line  used  in  the  analysis, 
Aq'   is  not  zero  and  positive  lover  surface  pressures  cannot  possibly  be 
realized  near  the  leading  edge  unless  the  angle  of  attack  is  increased  at 
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least  to  the  point  where  a  -  ^  =  0.   Since  Aj^  =  0   (with  n  >  2)  and 

A]_  sin  9  +  A2  sin  29  is  everywhere  positive  in  the  interval  0  S  9  ^  jt, 

the  condition  a  -  -^  =  0  is  sufficient  to  specify  positive  pressures 

o 
over  the  entire  chord  of  the  hydrofoil. 

A  convenient  way  of  treating  the  circular-arc  section  to  make  it 
comparahle  to  other  low-drag  sections  is  to  reorient  its  reference  line 

so  that  a,Q  =  ^,    that  is,   Aq'  =  0.   This  orientation  then  corresponds 
to  a  =  0.   At  this  design  angle  of  attack  la  =  0,  clq  =  ^j  the  hydro- 
dynamic  efficiency  Cj^JC-q     of  a  two-dimensional  circular-sire  section  as 
computed  from  linear  theory  is  ^  —2—,  which  is  almost  as  good  as  the 


low-drag  hydrofoil  selected  in  reference  1,  {  —^  =  ^ 


1^  2Cl      /Cl 

"»  j-^ -P  j-^  %^  r-i  v^  r»  *~i     I  '  ^^ 


VCd   16  2Cl 
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APPENDIX  B 

SAMPLE  CALCULATION  OF  LIFT,  DRAG,  AND  CENTER  OF  PRESSURE 
OF  A  CAMBERED  LIFTING  SURFACE  OPERATING  AT  FINITE 
DEPTH  AND  ZERO  CAVITATION  NUMBER 


A  sample  calculation  of  the  lift,  drag,  and  center  of  pressiire  is 
presented  for  a  Tiilln-Burkart  section  with  A-i_  =  0.2  fc^  d  ~  0-392), 
having  an  aspect  ratio  of  1  and  operating  at  an  angle  of  attack  of  12° 
and  depth  of  submersion  of  O.O7I  chord. 


LIFT  COEFFICIENT 


Step  1 

For  the  Tulin-Burkart  hydrofoil  section   (Ap  =   -  -ttI,   from  eq.ua- 
:ion   (7),  V''  2y 

o-coo  =  ^1-^  =  ^=  0-25 


Step  2 
From  figure  11  at  d/c  =  O.O7I, 


°^C,oo 


=  0.718 


Therefore, 


a„  =  (0.718) (0.25)  =  0.18 
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Step  5 

By  assuming  that     C^  1  =  0.25     and     t  =  0.12     ajid  by  using  equa- 
tion  (1+6),  ' 


Oi  =  1,12  ^^^  =  0.088 


and 


a  +  ttc  -  ai  =  -^=^  +   0.18  -  0.088  =  0.501  radian  =  I7.2 
^    ^   57.5 


From  figure  2  for  d/c  =  O.O7I  and  a  =  17.2°,  it  is  found  that  m  =  I.63; 
therefore,  from  reference  11,   t  =  0.12  as  assumed.   By  using  equa- 
tion (50), 

1  ,   ,  ,  ,     ,  cos  12° 

Cl,i  =  |(i.65) (0.501)       =  0.251 

'     ^  COS  -Ll.d. 

This  should  check  with  the  original  assumption  -  if  not,  repeat  step  5 
with  better  approximation  for  C^  and  t. 

Step  1+ 
For  the  Tulin-Burkart  section  given  by  equation  (12)  with  Ai  =  0.2, 


^(1)  =-^ 


Therefore, 


and 


ao  =  tan"^.035  =  1.92° 


a'  =  12°  +  1.92°  =  15.92^ 


a^  =  0.18  radian  =  10.3° 
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or 


a'  +  a^  =  2U.22° 


By  using  equation  (^8), 

Ct  „  =  i  0.88  sin2(2U.22)cos  15-92°  =  O.O72 
Thus,  the  required  lift  coefficient  is 

Cl  =  Cl  "L  +  Cl  (,  =  0.251  +  0.072  =  0.523 

DRAG  COEFFICIENT 

Step  1 
From  step  5  in  the  lift-coefficient  calculation, 

a  +  a^  ~  '^i  ~  17  •  2° 

For  d/c  =  0.071,  m  =  I.65;  and  for  d/c  =  00,  m  =  1.2.  Therefore,  from 
equation  (61) 

s   1.2  2 

From  figure  2  "by  using  the  a  =  0  line  with  m  =  2.15,  it  can  be  foimd 
that 

(d/c)g  =  0.66 
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Step  2 

If  a  (the  distance  to  airfoil  center  of  pressure)  is  assumed  to 
"be  0.57,  then  from  figure  9  for  a  =  O.37  and  d/c  =  0.66,  the  following 
values  are  found: 


Bo 

= 

0. 

56 

Bl 

= 

0. 

.15 

B2 

^w 

0 

From  equation  (35)  for  Aq  =  a  -  a^^  =  0.121,   A-[_  =  0.2  and  k^  =   -0.1. 
The  values  of  Cq,  C-^,    and  C2  may  be  determined  as 


Cn  =  0.212 


c-,_  =  0.181 

Cg  =  -0.1 


and  from  equation  (29),   a  =  O.567,  which  checks  with  the  assumed  value. 
If  the  resulting  value  of  a  differs  enough  from  the  assumed  value  to 
affect  the  values  of  the  B  coefficients,  step  2  should  be  repeated. 


Step  3 
From  equation  (60)  the  value  of  e   is  determined  as 

e  =  0.0013  radian  =  0.065° 

Step  k 
From  equation   (5l)> 

Cl  1  =  0.251  Cl  c  =  0-072  Cf  =  O.OOU 

Therefore, 

Cd  =  0.251  tan(l2.065°)    +  O.O72  tan(l3.92°)    +  O.OOli  =  O.O75 
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CENTER  OF  PRESSURE 

Step  1 

From  step  2  of  the  drag-coefficient  calculation, 

Cq  =  0.212     Ci  =  0.181     C2  =  -0.1 
Thus,  from  equation  (62), 

^c.p.,l=  0-^5c 

Step  2 

From  equation  (6h)    and  "by  using 

%1  -  0-251 
and 

%c  =  0-072 

it  is  found  that 


^c.p.   =  0.i+5c 
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(a)  Supercavltating  flow  at  great  depth  of  submersion. 


(b)  Supercavltating  or  ventilated  flow  near  the  free  surface. 
Figure  1.-  Definition  sketch. 
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(a)  Hydrofoil,  Z-plane. 


Transformed  upper 
cavity  surface 


Trans  formed 
water  surface 

icy  =  d/2 


Center  of  pressure 


n.O s.  J 


Transformed  lower 
cavity  surface 


("b)  Equivalent  airfoil,  Z-plane. 


Figure  5--  Linearized  boundary  conditions  in  hydrofoil  and  equivalent 

airfoil  planes . 
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Equivalent  airfoil ,   S  — plane 

Figure  6.-  Linearized  model  for  calculating  effect  of  depth  of  submer- 
sion on  lift  coefficient  of  flat  plate  (a  =  O). 
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Figiire  7'-  Comparison  of  linearized  solution  with  exact  soultion  for 
effect  of  depth  of  submersion  on  lift  coefficient  of  a  flat  plate 
(a  =  0). 
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Water  surface 


Stagnation  streamline 


d  =  5  -  %  =  05 


(a)   6/c  =  00. 


Free  surface 


(b)   5/c  =  0. 


(c)    0  <  S/c  <  CO. 

Figure  12.-  Location  and  direction  of  stagnation  line  on  flat  plate 
normal  to  direction  of  motion. 
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face 


(a)  Hydrofoil,  Z-plane. 


Transformed  upper 
cavity  surface 


(b)  Airfoil,   Z  =  -  \fz'  plane. 


Figure  I5.-  Hydrofoil  and  equivalent  airfoil  at  large  angles  of  attack 

and  small  depths  of  submersion. 
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L- 57- 27^5 
Figure  I5.-  Test  set-up  showing  cambered  hydrofoil  with  aspect  ratio 
of  1  and  balance  attached  to  towing  carriage. 
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Figure  16.-  System  of  axes. 
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(b)    Front  elevation. 
Figure   17--   Schematic  drawing  of  spray-thickness   gauge. 
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V  =  20  fps 


V  =  kO   fps 


V  =  60  fps  V  =  80  fps 

(a)   a  =  8°. 


V  =  10  fps 


V  =  30  fps 


(b)   a  =  10°. 


V  =  20  fps 


V  =  1+0  fps 


L- 57- 27^0 


Figure  18.- 


Flow  about  the  flat  lifting  surface,  depth  of  submersion  is 

0.5  inch. 
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V  =  10  fps 


V  =  30  fps 


V  =  20  fps 


V  =  kO   fps 


(c)   a  =  8°.   Air  supplied  at  0.012  pound  per  second  for  V  =  10,  20, 

30  feet  per  second. 


Figure  I8.-  Concluded. 


V  =  20  fps 


V  =  60  fps 


V  =  UO  fps 


V  =  80  fps 


(a)   a  =  10°. 


L- 57- 27^4-1 


Figure  19--  Flow  about  cambered  lifting  surface,  depth  of  submersion  is 

0.5  inch. 
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V  =  20  f  ps 


V  =  i+0  fps 


V  =  59  fps  V  =  80  fps 

(b)   a  =  12°. 


V  =  10  fps 


V  =  20  fps 


V  =  30  fps  V  =  40  fps 

(c)   a  =  16°.  L-57-2Ti^2 

Figiire  I9.-  Concluded. 
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(a)  Lift. 
Figure  22.-  Characteristics  of  ventilated  flat  lifting  surface. 
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Figure  22.-  Continued. 
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(c)  Pitching  moment  about  leading  edge. 
Figure  22.-  Concluded. 
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(a)  Lift. 
Figure  23.-  Characteristics  of  ventilated  cambered  lifting  surface. 


NACA  EM  L57II6 


79 


Depth, 
in. 

O  1/2 
D  1/2 
O  1/2 
A  1/2 
t^  1/2 
b  1/2 
Q  1/2 

c:^  0 
a  0 
-^   0 


Of, 
deg 

20 
18 
16 
114 
12 
10 

g 
10 

& 

6 


Venti lation 

Vortex 

Vortex 

Vortex 

Vortex 

Vortex  _ 

Vortex,  wire  on  L.E. 

Forced  air,  wire  on  L.E. 

Vortex 

Vortex 


10 


20 


30      MO      50 
Speed,  fps 


70 


90 


(b)  Drag. 
Figure  25.-  Continued. 
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Figure  23.-  Concluded. 
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Figure  24.-  Characteristics   of  ventilated  flat  plate  at     a  =  l6°  and  20° 
for  depths   of  submersion  of  0,    1.0,    I.5,   and  2.0  inches. 
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Figure  2h.-   Continued. 
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Figure  33--  Effect  of  cavitation  number  on  lift  coefficients  of  flat 
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aspect  ratio  of  1. 
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Figure  3U.-  Comparison  of  air-  and  vapor -filled  cavities  for  depth  of 
submersion  of  6  inches  and  a  of  20°. 
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